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ABSTRACT 

We present spectral and photometric observations of 10 type la supernovae (SNe la) in the 
redshift range 0.16 < z < 0.62. The luminosity distances of these objects are determined by 
methods that employ relations between SN la luminosity and light curve shape. Combined with 
previous data from our High-Z Supernova Search Team (Garnavich et al. 1998; Schmidt et al. 
1998) and Riess et al. (1998a), this expanded set of 16 high-redshift supernovae and a set of 
34 nearby supernovae are used to place constraints on the following cosmological parameters: 
the Hubble constant (Hq), the mass density (fijvf), the cosmological constant (i.e., the vacuum 
energy density, Oa), the deceleration parameter (go), and the dynamical age of the Universe (to). 
The distances of the high-redshift SNe la are, on average, 10% to 15% farther than expected 
in a low mass density (CIm = 0.2) Universe without a cosmological constant. Different light 
curve fitting methods, SN la subsamples, and prior constraints unanimously favor eternally 
expanding models with positive cosmological constant (i.e., Q,\ > 0) and a current acceleration 
of the expansion (i.e., qo < 0). With no prior constraint on mass density other than D,m > 0, 
the spectroscopically confirmed SNe la are statistically consistent with qo < at the 2.8er 
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and 3.9<7 confidence levels, and with Ct\ > at the 3. Oct and 4. Oct confidence levels, for two 
different fitting methods respectively. Fixing a "minimal" mass density, J7m = 0.2, results in 
the weakest detection, > at the 3. Oct confidence level from one of the two methods. For 
a flat-Universe prior (Q M + Cl\ = 1), the spectroscopically confirmed SNe la require Sl\ > 
at 7ct and 9ct formal significance for the two different fitting methods. A Universe closed by 
ordinary matter (i.e., Om = 1) is formally ruled out at the 7ct to 8ct confidence level for the 
two different fitting methods. We estimate the dynamical age of the Universe to be 14.2 ±1.5 
Gyr including systematic uncertainties in the current Cepheid distance scale. We estimate the 
likely effect of several sources of systematic error, including progenitor and metallicity evolution, 
extinction, sample selection bias, local perturbations in the expansion rate, gravitational lensing, 
and sample contamination. Presently, none of these effects reconciles the data with fi A = and 
qo > 0. 

subject headings: supernovae: general ; cosmology:observations 

1. Introduction 

This paper reports observations of 10 new high-redshift type la supernovae (SNe la) and the values of 
the cosmological parameters derived from them. Together with the four high-redshift supernovae previously 
reported by our High-Z Supernova Search Team (Schmidt et al. 1998; Garnavich et al. 1998) and two 
others (Riess et al. 1998a), the sample of 16 is now large enough to yield interesting cosmological results of 
high statistical significance. Confidence in these results depends not on increasing the sample size but on 
improving our understanding of systematic uncertainties. 

The time evolution of the cosmic scale factor depends on the composition of mass-energy in the 
Universe. While the Universe is known to contain a significant amount of ordinary matter, Qm, which 
decelerates the expansion, its dynamics may also be significantly affected by more exotic forms of energy. 
Pre-eminent among these is a possible energy of the vacuum Einstein's "cosmological constant," 

whose negative pressure would do work to accelerate the expansion (Carroll, Press, & Turner 1992; Schmidt 
et al. 1998). Measurements of the redshift and apparent brightness of SN la of known intrinsic brightness 
can constrain these cosmological parameters. 

1.1. The High-Z Program 

Measurement of the elusive cosmic parameters Hm and Oa through the redshift-distance relation 
depends on comparing the apparent magnitudes of low-redshift SNe la with those of their high-redshift 
cousins. This requires great care to assure uniform treatment of both the nearby and distant samples. 

The High-Z Supernova Search Team has embarked on a program to measure supernovae at high 
redshift and to develop the comprehensive understanding of their properties required for their reliable use 
in cosmological work. Our team pioneered the use of supernova light curve shapes to reduce the scatter 
about the Hubble line from a k, 0.4 mag to a « 0.15 mag (Hamuy et al. 1996a, b, 1995; Riess, Press 
& Kirshner 1995, 1996a). This dramatic improvement in the precision of SNe la as distance indicators 
increases the power of statistical inference for each object by an order of magnitude and sharply reduces 
their susceptibility to selection bias. Our team has also pioneered methods for using multi-color observations 
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to estimate the reddening to each individual supernova, near and far, with the aim of minimizing the 
confusion between effects of cosmology and dust (Riess, Press, & Kirshner 1996a; Phillips et al. 1998). 
Because the remaining scatter about the Hubble line is so small, the discussion of the Hubble constant from 
low-rcdshift SNe la has already passed into a discussion of the best use of Cepheid distances to galaxies that 
have hosted SNe la (Saha et al. 1997; Kochanek 1997; Madore & Freedman 1998; Riess, Press, & Kirshner 
1996a; Hamuy et al. 1996a; Branch 1998). As the use of SNe la for measuring Qm and I^a progresses from 
its infancy into childhood, we can expect a similar shift in the discussion from results limited principally by 
statistical errors to those limited by our depth of understanding of SNe la. 

Published high-redshift SN la data is a small fraction of the data in hand both for our team and for the 
Supernova Cosmology Project (Perlmuttcr et al. 1995, 1997, 1998). Now is an opportune time to spell out 
details of the analysis, since further increasing the sample size without scrupulous attention to photometric 
calibration, uniform treatment of nearby and distant samples, and an effective way to deal with reddening 
will not be profitable. Besides presenting results for four high-z supernovae, we have published details of 
our photometric system (Schmidt et al. 1998) and stated precisely how we used ground-based photometry 
to calibrate our Hubble Space Telescope (HST) light curves (Garnavich et al. 1998). In this paper, we spell 
out details of newly-observed light curves for 10 objects, explain the recalibration of the relation of light 
curve shape and luminosity for a large low-redshift sample, and combine all the data from our team's work 
to constrain cosmological parameters. We also evaluate how systematic effects could alter the conclusions. 
While some comparison with the stated results of the Supernova Cosmology Project (Perlmutter et al. 
1995, 1997, 1998) is possible, an informed combination of the data will have to await a similarly detailed 
description of their measurements. 

1.2. A Brief History of Supernova Cosmology 

While this paper emphasizes new data and constraints for cosmology, a brief summary of the subject 
may help readers connect work on supernovae with other approaches to measuring cosmological parameters. 

Empirical evidence for SNe I presented by Kowal (1968) showed that these events had a well-defined 
Hubble diagram whose intercept could provide a good measurement of the Hubble constant. Subsequent 
evidence showed that the original spectroscopic class of Type I should be split (Doggett & Branch 1985; 
Uomoto & Kirshner 1985; Wheeler & Lcvreault 1985; Wheeler & Harkness 1986; Porter & Filippenko 1987). 
The remainder of the original group, now called Type la, had peak brightness dispersions of 0.4 mag to 0.6 
mag (Tammann & Lcibundgut 1990; Branch & Miller 1993; Miller & Branch 1990; Delia Vallc & Panagia 
1992; Rood 1994; Sandage & Tammann 1993; Sandage et al. 1994). Theoretical models suggested that 
these "standard candles" arose from the thermonuclear explosion of a carbon-oxygen white dwarf that had 
grown to the Chandrasckhar mass (Hoyle & Fowler 1960; Arnett 1969; Colgate & McKee 1969). Because 
SNe la are so luminous {Mb ~ —19.5 mag), Colgate (1979) suggested that observations of SNe la at z ss 1 
with the forthcoming Space Telescope could measure the deceleration parameter, q . 

From a methodical CCD-based supernova search that spaced observations across a lunation and 
employed prescient use of image-subtraction techniques to reveal new objects, Hansen, N0rgaard-Nielsen, & 
Jorgensen (1987) detected SN 1988U, a SN la at z = 0.31 (N0rgaard-Nielsen et al. 1989). At this redshift 
and distance precision (a « 0.4 to 0.6 mag), ~ 100 SNe la would have been needed to distinguish between 
an open and closed Universe. Since the Danish group had already spent two years to find one object, it was 
clear that larger detectors and faster telescopes needed to be applied to this problem. 
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Evidence of systematic problems also lurked in supernova photometry so that merely increasing the 
sample would not be adequate. Attempts to correct supernova magnitudes for reddening by dust (Branch 
& Tammann 1992) based on the plausible (but incorrect) assumption that all SNe la had the same intrinsic 
color had the unfortunate effect of increasing the scatter about the Hubble line or alternately attributing 
bizarre properties to the dust absorbing SN la light in other galaxies. In addition, well-observed supernovae 
such as SN 1986G (Phillips et al. 1987; Cristiani ct al. 1992), SN 1991T (Filippenko et al. 1992a; Phillips 
ct al. 1992; Ruiz-Lapuente et al. 1992), and SN 1991bg (Filippenko et al. 1992b; Leibundgut et al. 1993; 
Turatto et al. 1996) indicated that large and real inhomogeneity was buried in the scatter about the Hubble 
line. 

Deeper understanding of low-redshift supernovae greatly improved their cosmological utility. Phillips 
(1993) reported that the observed peak luminosity of SNe la varied by a factor of 3. But he also showed 
that the decrease in B brightness in the 15 days after peak (Ami 5 (5)) was a good predictor of the SN la 
luminosity with slowly declining supernovae more luminous than those that fade rapidly. 

A more extensive database of carefully and uniformly observed SNe la was needed to refine the 
understanding of SN la light curves. The Calan/Tololo survey (Hamuy et al. 1993a) made a systematic 
photographic search for supernovae between cycles of the full moon. This search was extensive enough to 
guarantee the need for scheduled follow-up observations, which were supplemented by the cooperation of 
visiting observers, to collect well-sampled light curves. Analysis of the Calan/Tololo results generated a 
broad understanding of SNe la and demonstrated their remarkable distance precision (after template fitting) 
of d ~ 0.15 mag (Hamuy et al. 1995, 1996a, b,c,d). A parallel effort employed data from the Calan/Tololo 
survey and from the Harvard-Smithsonian Center for Astrophysics (CfA) to develop detailed empirical 
models of SN la light curves (Riess, Press, & Kirshner 1995; Riess 1996). This work was extended into the 
Multi-Color Light Curve Shape (MLCS) method which employs up to 4 colors of SN la photometry to yield 
excellent distance precision (w 0.15 mag) and a statistically valid estimate of the uncertainty for each object 
with a measurement of the reddening by dust for each event (Riess, Press, & Kirshner 1996a; Appendix of 
this paper). This work has also placed useful constraints on the nature of dust in other galaxies (Riess, 
Press, & Kirshner 1996b). 

The complete sample of nearby SNe la light curves from the Calan/Tololo and CfA samples provides 
a solid foundation from which to extend the redshift-distancc relation to explore cosmological parameters. 
The low-redshift sample used here has 34 SNe la with z < 0.15. 

Since the high-redshift observations reported here consumed large amounts of observing time at the 
world's finest telescopes, we have a strong incentive to find efficient ways to use the minimum set of 
observations to derive the distance to each supernova. A recent exploration of this by Riess et al. (1998a) 
is the "Snapshot" method which uses only a single spectrum and a single set of photometric measurements 
to infer the luminosity distance to a SN la with ~ 10% precision. In this paper, we employ the snapshot 
method for six SNe la with sparse data, but a shrewdly designed program that was intended to use the 
snapshot approach could be even more effective in extracting useful results from slim slices of observing 
time. 

Application of large-format CCDs and sophisticated image analysis techniques by the Supernova 
Cosmology Project (Perlmutter et al. 1995) led to the discovery of SN 1992bi (z — 0.46) followed by 6 more 
SNe la at z « 0.4 (Perlmutter et al. 1997). Employing a correction for the luminosity/light-curve shape 
relation (but none for host galaxy extinction), comparison of these SNe la to the Calan/Tololo sample gave 
an initial indication of "low" Q\ and "high" Qm- = 0. 061^34 for a flat Universe and f^M = 0.88^;^ 
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for a Universe without a cosmological constant = 0). The addition of one very high-rcdshift (z = 0.83) 
SN la observed with HST had a significant effect on the results: Qa = 0.4 ± 0.2 for a flat Universe, and 
M = 0.2 ± 0.4 for a Universe with = 0. (Perlmutter et al. 1998). This illustrates how young and 
volatile the subject is at present. 

1.3. This Paper 

Our own High-Z Supernova Search Team has been assiduously discovering high-redshift supernovae, 
obtaining their spectra, and measuring their light curves since 1995 (Schmidt et al. 1998). The goal is to 
provide an independent set of measurements that uses our own techniques and compares our data at high 
and low rcdshifts to constrain the cosmological parameters. Early results from 4 SNc la (3 observed with 
HST) hinted at a non-negligible cosmological constant and "low" CIm, but were limited by statistical errors: 
Ha = 0.65 ± 0.3 for a flat Universe, VLm = —0.1 ± 0.5 when Oa = (Garnavich et al. 1998). Our aim in 
this paper is to move the discussion forward by increasing the data set from four high-redshift SNe to 16, to 
spell out exactly how we have made the measurement, and to consider various possible systematic effects. 

In §2 we describe the observations of the SNe la including their discovery, spectral identification, 
photometric calibration, and light curves. We determine the luminosity distances (including X-corrections) 
via two methods, MLCS and a template fitting method (Ami^(B)) , as explained in §3. Statistical inference 
of the cosmological parameters including H , Q M , f^, q a , t a , and the fate of the Universe is contained in 
§4. Section 5 presents a quantitative discussion of systematic uncertainties which could affect our results: 
evolution, absorption, selection bias, a local void, weak lensing, and sample contamination. Our conclusions 
are summarized in §6. 



2. Observations 
2.1. Discovery 

We have designed a search program to find supernovae in the redshift range 0.3 < z < 0.6 with the 
purpose of measuring luminosity distances to constrain cosmological parameters (Schmidt et al. 1998). 
Distances with the highest precision are measured from SNe la observed before maximum brightness and 
in the redshift range of 0.35 < z < 0.55, where our set of custom passbands measures the supernova light 
emitted in rest- frame B and V. By imaging fields near the end of a dark run, and then again at the 
beginning of the next dark run, we ensure that the newly discovered supernovae are young (N0rgaard-Nielsen 
et al. 1989; Hamuy et al. 1993a; Perlmutter et al. 1995). Observing a large area and achieving a limiting 
magnitude of m R w 23 mag yields many SN la candidates in the desired redshift range (Schmidt et al. 
1998). By obtaining spectra of these candidates with 4-m to 10-m telescopes, we can identify the SNe la 
and confirm their youth using the spectral feature aging technique of Riess et al. (1997). 

The 10 new SNc la presented in this paper (SN 1995ao, SN 1995ap, SN 1996E, SN 1996H, SN 19961, 
SN 1996J, SN 1996K, SN 1996R, SN 1996T, and SN 1996U) were discovered using the CTIO Blanco 4 m 
telescope with the facility prime-focus CCD camera as part of a 3-night program in 1995 Oct-Nov and a 
6-night program in 1996 Feb-Mar. This instrument has a pixel scale of 0.43", and the TEK 2048x2048 pixel 
CCD frame covers 0.06 square degrees. In each of the search programs, multiple images were combined 
after removing cosmic rays, differenced with "template" images, and searched for new objects using the 
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prescription of Schmidt et al. (1998). The data on 1995 Oct 27 and 1995 Nov 17 were gathered under 
mediocre conditions with most images having seeing worse than 1.5". The resulting differenced images 
were sufficient to find new objects brighter than = 22.5 mag. The data acquired in 1996 had better 
image quality (~ 1.5"), and the differenced images were sufficient to uncover new objects brighter than 
run = 23 mag. 

In total, 19 objects were identified as possible supernovae — 2 new objects were detected on both 1995 
November 17 and 1995 November 29, 5 new objects on 1996 Feb 14-15, 2 on 1996 Feb 20-21, and 8 on 1996 
Mar 15-16 (Kirshner et al. 1995; Garnavich et al. 1996a, b). 

2.2. Data 

Spectra of the supernova candidates were obtained to classify the SNe and obtain redshifts of their 
host galaxies. For this purpose, the Keck telescope, Multiple-Mirror Telescope (MMT), and the European 
Southern Observatory 3.6-m (ESO 3.6-m) were utilized following the Fall 1995 and Spring 1996 search 
campaigns. Some galaxy redshifts were obtained with the Keck telescope in the Spring of 1998. 

The Keck spectra were taken with the Low Resolution Imaging Spectrograph (LRIS; Oke et al. 1995), 
providing a resolution of 6 A full width at half maximum (FWHM). Exposure times were between 3x900 
seconds and 5x900 seconds, depending on the candidate brightness. 

The MMT spectra were obtained with the Blue Channel spectrograph and 500 lines/mm grating giving 
a resolution of 3.5 A FWHM. Exposure times were 1200 seconds and repeated five to seven times. The 
MMT targets were placed on the slit using an offset from a nearby bright star. 

The ESO 3.6-m data were collected with the ESO Faint Object Spectrograph Camera (EFOSC1) at a 
nominal resolution of 18 A FWHM. Single 2700 second exposures were made of each target. 

Using standard reduction packages in IRAF, the CCD images were bias subtracted and divided by a 
flat-field frame created from a continuum lamp exposure. Multiple images of the same object were shifted 
where necessary and combined using a median algorithm to remove cosmic ray events. For single exposures, 
cosmic rays were removed by hand using the IRAF/IMEDIT routine. Sky emission lines were problematic, 
especially longward of 8000 A . The spectra were averaged perpendicular to the dispersion direction, and 
that average was subtracted from each line along the dispersion. However, residual noise from the sky lines 
remains. The one-dimensional spectra were then extracted using the IRAF/APSUM routine and wavelength 
calibrated either from a comparison lamp exposure or the sky emission lines. The flux was calibrated using 
observations of standard stars and the IRAF/ONEDSTDS database. 

The candidates were classified from visual inspection of their spectra and comparison with the spectra 
of well-observed supernovae (see §5.7). In all, 10 of the candidates were SNe la, 1 was a SN II, and 2 were 
active galactic nuclei or SNe II (Kirshner et al. 1995; Garnavich et al. 1996a, b). The remaining 6 candidates 
were observed, but the spectra did not have sufficient signal to allow an unambiguous classification. 
The identification spectra for the 10 new SNe la are summarized in Table 1 and shown in Figure 1. In 
addition we include the spectral data for 3 previously analyzed SNe: SN 1997ce, SN 1997cj, and SN 1997ck 
(Garnavich et al. 1998). The spectral data for SN 1995K are given by Schmidt et al. (1998). The spectrum 
of SN 1997ck shows only a [O II] emission line at 7328.9 A in four separate exposures (Garnavich et al. 
1998). The equivalent R band magnitude of the exposure was 26.5 which is more than 1.5 magnitudes 
dimmer than the supernova would have been in R, suggesting that the SN was not in the slit when the host 
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galaxy was observed. 

Most of the host galaxies showed emission lines of [O II], [O III], or Ha in the spectrum, and the 
rcdshift was easily measured for these. For the remainder, the redshift was found by matching the broad 
features in the high-redshift supernovae to those in local supernova spectra. The intrinsic dispersion in the 
expansion velocities of SNe la (Branch et al. 1988; Branch & van den Bergh 1993) limits the precision of 
this method to la ~ 2500 km s _1 independent of the signal-to-noise ratio of the SN spectrum. The method 
used to determine the redshift for each SN is given in Table 1. 

Following the discovery and identification of the SNe la, photometry of these objects was obtained from 
observatories scheduled around the world. The SNe were primarily observed through custom passbands 
designed to match the wavelength range closest to rest- frame Johnson B and V passbands. Our "£?45," 
"F45," ".B35," and "F35" filters are specifically designed to match Johnson B and V redshifted by z = 0.45 
and z — 0.35, respectively. The characteristics of these filters are described by Schmidt et al. (1998). A few 
observations were obtained through standard bandpasses as noted in Table 2 where we list the photometric 
observations for each SN la. 

Photometry of local standard stars in the supernova fields in the £>35, F35, £?45, F45 (or "supernova") 
photometric system were derived from data taken on three photometric nights. The method has been 
described in Schmidt et al. (1998) but we summarize it here. The supernova photometric system has 
been defined by integrating the fluxes of spectrophotometric standards from Hamuy et al. (1994) through 
the supernova bandpass response functions (based on the filter transmissions and a typical CCD quantum 
efficiency function) and solving for the photometric coefficients which would yield zero color for these stars 
and monochromatic magnitudes of 0.03 for Vega. 

This theoretically defined photometric system also provides transformations between the Johnson/Kron- 
Cousins system and the supernova system. We use theoretically derived transformations to convert the 
known V, R, and I magnitudes of Landolt (1992) standard fields into -B35, V35, 545, V45 photometry. 

On nights which are photometric, we observe Landolt standard fields with the B35, V35, BA5, V45 
filters and measure the stars' instrumental magnitudes from apertures large enough to collect all the 
stellar light. We then derive the transformation from the supernova system to the instrumental system 
as a function of the instrumental magnitudes, supernova system colors, and observed airmass. Because 
our theoretical response functions are very similar to the instrumental response functions, our measured 
color coefficients were small, typically less than 0.02 mag per mag of £?45 — VA5 or B35 — 1^35. These 
long wavelength filters also reduced the effect of atmospheric extinction (compared to B and V). Typical 
extinction coefficients were 0.11, 0.09, 0.07, and 0.06 mag per airmass for 535, £>45, F35, and V45, 
respectively. 

Isolated stars on each supernova frame were selected as local standards. The magnitudes of the local 
standards were determined from the transformation of their instrumental magnitudes, measured from 
similarly large apertures. The final transformed magnitudes of these local standards, averaged over three 
photometric nights, is given in Table 3. The locations of the local standards and the SNe are shown in 
Figure 2. The uncertainties in the local standards' magnitudes are the quadrature sum of the uncertainty 
(dispersion) of the instrumental transformations (typically 0.02 mag) and the individual uncertainties from 
photon (Poisson) statistics. The dispersion in the instrumental transformation quantifies the errors due to 
imperfect flat-fielding, small changes in the atmospheric transparency, incomplete empirical modeling of 
the response function, and seeing variations. This uncertainty is valid for any single observation of the local 
standards. 
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To measure the brightness of the supernovae free from host galaxy contamination, we obtained deep 
images of the hosts a year after, or months before, the discovery of the SNe. These images were used to 
subtract digitally a host's light from the supernova's light, leaving only the stellar point spread function 
(PSF). The algorithms employed to match the resolution, intensity, and coordinate frames of images prior 
to subtraction are described in Schmidt et al. (1998). The brightness of the SNe in these uncrowded fields 
was then measured relative to the calibrated local standard stars in the field by fitting a model of a PSF 
to the stars and supernova using the DoPHOT algorithm (Schmidt et al. 1998; Mateo & Schechter 1989; 
Schechter, Mateo, & Saha 1993). 

Systematic and statistical components of error were evaluated by measuring the brightness of artificial 
stars added to the subtracted frames. These artificial stars had the same brightness and background as the 
measured SNe (Schmidt et al. 1998). The "systematic" error was measured from the difference in the mean 
magnitude of the artificial stars before and after the image processing (i.e., alignment, scaling, "blurring," 
and subtracting). The systematic errors were always < 0.1 mag and were of either sign. Any significant 
systematic error is likely the result of a mismatch in the global properties of the template image and SN 
image based on only examining a local region of the two images. A correction based on the systematic 
error determined from the artificial stars was applied to the measured SN magnitude to yield an unbiased 
estimate of the SN magnitude. The dispersion of the recovered artificial magnitudes about their mean was 
assigned to the statistical uncertainty of the SN magnitude. 

The supernova PSF magnitudes were transformed to the -B35, V35, BA5, VA5 system using the local 
standard magnitudes and the color coefficients derived from observations of the Landolt standards. The 
final SN light curves are the average of the results derived from 5 or 6 local standards, weighted by the 
uncertainty of each local standard star. The light curves are listed in Table 4 and displayed in Figure 3. 
The SN magnitude errors are derived from the artificial star measurements as described above. 

The small color and atmospheric extinction coefficients give us confidence that the supernova 
photometry accurately transformed to the B35, V35, £?45, V45 system. However, it is well known that a 
nonstellar flux distribution can produce substantial systematic errors in supernova photometry (Menzies 
1989). We have anticipated this problem by using identical filter sets at the various observatories and by 
defining our photometric system with actual instrumental response functions. To measure the size of this 
effect on our SN photometry, we have calculated the systematic error incurred from the differences in the 
instrumental response functions of different observatories we employed. Spectrophotometric calculations 
from SN la spectra using various instrumental response functions show that the expected differences are 
less than 0.01 mag and can safely be ignored. 

3. Analysis 
3.1. i^-corrections 

A strong empirical understanding of SN la light curves has been garnered from intensive monitoring of 
SNe la at z < 0.1 through B and V passbands (Hamuy et al. 1996b; Riess 1996; Riess et al. 1998b; Ford 
et al. 1993; Branch 1998, and references therein). We use this understanding to compare the light curves 
of the high-redshift and low-redshift samples at the same rest wavelength. By a judicious choice of filters, 
we minimize the differences between B and V rest-frame light observed for distant SNe and their nearby 
counterparts. Nevertheless, the range of rcdshifts involved makes it difficult to eliminate all such differences. 
We therefore employ "if-corrections" to convert the observed magnitudes to rest-frame B and V (Oke & 
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Sandage 1968; Hamuy et al. 1993b; Kim, Goobar, & Perlmutter 1996; Schmidt et al. 1998). 

The cross-band if-correction for SNe la has been described as a function of the observed and rest-frame 
filter transmissions, the redshift of the supernova, and the age of the supernova (see equation 1 of Kim, 
Goobar, & Perlmutter 1996). Such a if -correction assumes that the spectral energy distribution of all 
SNe la of a given age is homogeneous, yet it has been shown (Pskovskii 1984; Phillips et al. 1987, 1993; 
Leibundgut et al. 1993; Nugent et al. 1995; Riess, Press, & Kirshner 1996a; Phillips et al. 1998; Lira 
1995; Appendix of this paper) that at a given age, the colors of SNe la exhibit real variation related to the 
absolute magnitude of the supernova. 

A variation in SN la color, at a fixed phase, could have dire consequences for determining accurate 
if -corrections. An appropriate if-correction quantifies the difference between the supernova light which 
falls into a standard passband (e.g., B) at z = and that which falls into the filters we employ to observe 
a redshifted SN la. Differences in SN la color, at a fixed phase, would alter the appropriate if-correction. 
We need to know the color of each supernova to determine its if-correction precisely. Differences in SN la 
color can arise from interstellar extinction or intrinsic properties of the supernova such as a variation in 
photospheric temperature (Nugent et al. 1995). 

Nugent et al. (1998) have shown that, to within 0.01 mag, both the effects of extinction and intrinsic 
variations on the SN la spectral energy distribution near rest-frame B and V and hence on the if-correction 
can be reproduced by application of a Galactic reddening law (Cardclli, Clayton, & Mathis 1989) to the 
spectra. The difference in color, at a given age, between an individual SN la and a fiducial SN la is 
quantified by a color excess, Eb-v, and determines the effects of either extinction or intrinsic variation on 
the spectra and observed colors of the SNe. For most epochs, filter combinations, and redshifts, the variation 
of the if-correction with the observed variations of color excess is only 0.01 to 0.05 mag. For redshifts 
which poorly match the rest-frame wavelengths to the observed wavelengths, the custom if -correction for 
very red or very blue SNe la can differ from the standard if -correction by 0.1 to 0.2 mag. 

This prescription requires the age and observed color for each observation to be known before its 
if-correction can be calculated. The age is best determined from fitting the light curve's time of maximum. 
Yet we must use the if-correction to determine the time of maximum and the true color of each epoch. 
This conundrum can be solved by iteratively converging to a solution by repeated cycles of if-correcting 
and empirical fitting of the light curves. Table 4 lists the final cross-band if -corrections we used to convert 
the observations to the rest-frame passbands. 

3.2. Light Curve Fitting 

As described in §1, empirical models for SNe la light curves which employ the observed correlation 
between light curve shape, luminosity, and color have led to significant improvements in the precision of 
distance estimates derived from SNe la (Hamuy et al. 1995, 1996a; Riess, Press, & Kirshner 1995, 1996a). 
Here we employ the MLCS method prescribed by Riess, Press, & Kirshner 1996a as reanalyzed in the 
Appendix and the template fitting method of Hamuy et al. (1995, 1996d) to fit the light curves in Table 4. 

The growing sample of well-observed SN la light curves (Hamuy et al. 1996b; Riess 1996; Riess et al. 
1998b; Ford et al. 1993) justifies refinements in the MLCS method that are described in the Appendix. 
These include a new derivation of the relation between light curve shape, luminosity, and color from SNe 
la in the Hubble flow using redshift as the distance indicator. In addition, this empirical description has 
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been extended to a second order (i.e., quadratic) relation between SN la luminosity and light curve shape. 
A more realistic a priori probability distribution for extinction has been utilized from the calculations of 
Hatano, Branch, & Deaton (1997). Further, we now quantify the residual correlations between observations 
of dissimilar time, passband, or both. The empirical model for a SN la light and color curve is still described 
by four parameters: a date of maximum (t), a luminosity difference (A), an apparent distance (/xs), and an 
extinction (As)- Due to the redshifts of the SN host galaxies we first correct the supernova light curves for 
Galactic extinction (Burstein & Heiles 1982), then determine host galaxy extinction. 

To treat the high and low redshift SNe la consistently, we restricted the MLCS fits to the nearby 
SNe la observations in B and V within 40 days after maximum brightness in the restframe. This is the 
age by which all high-rcdshift light curve observations ended. Because of this restriction, we also limited 
our consideration of nearby SNe la to those with light curves which began no later than ~ 5 days after B 
maximum. Although more precise distance estimates could be obtained for the nearby sample by including 
later data and additional colors, the nearby sample is large enough to determine the nearby expansion rate 
to sufficient precision. The parameters of the MLCS fits to 27 SNe la in the nearby Hubble flow (0.01 < z 
< 0.13; Hamuy ct al. 1996b; Riess et al. 1998b) arc given in Tabic 10. 

In Table 5 we list the parameters of the MLCS fits to six SN la light curves presented here (SNe 1996E, 
1996H, 19961, 1996J, 1996K, 1996U) and for three SNe la from our previous work (SNe 1995K, 1997ce, 
1997cj; Garnavich et al. 1998; Schmidt et al. 1998). We have placed all MLCS distances on the Cephcid 
distance scale using Cepheid distances to galaxies hosting photoelectrically observed SNe la: SN 1981B, 
SN 1990N, and SN 1972E (Saha et al. 1994, 1997; Riess, Press, & Kirshner 1996a). However, conclusions 
about the values of the cosmological parameters Qm, ^a, and go are independent of the distance scale. 

An additional supernova, SN 1997ck, was studied by Garnavich et al. (1998) in a galaxy with z=0.97. 
Its rest-frame B light curve was measured with the HST (see Figure 3). Although this object lacks a 
spectroscopic classification and useful color information, its light curve shape and peak luminosity are 
consistent with those of a typical SN la. Due to the uncertainty in this object's extinction and classification, 
we will analyze the SNe la distances both with and without this most distant object. 

We have also determined the distances to the same 27 nearby SNe la and the ten well-observed 
high-redshift events using a template fitting approach (Hamuy et al. 1995, 1996d). The maximum light 
magnitudes and the initial decline rate parameter Amis(S) for a given SN la are derived by comparing 
the goodness-of-fits of the photometric data to a set of six template SN la light curves selected to cover 
the full range of observed decline rates. The intrinsic luminosity of the SN is then corrected to a standard 
value of the decline rate (Ami5(B)=l.l) using a linear relation between Ami5(B) and the luminosities for 
a set of SNe la in the Hubble flow (Phillips et al. 1998). An extinction correction has been applied to these 
distances based on the measured color excess at maximum light using the relation between Ami5(B) and 
the unreddened SN la color at maximum light (Phillips et al. 1998). These extinction measurements employ 
the same Bayesian filter (in the Appendix) used for the MLCS fits. The final distance moduli are also on 
the Cepheid distance scale as described by Hamuy et al. (1996a) and Phillips et al. (1998). Parameters of 
these fits to the nearby and high-redshift SNe la are provided in Table 10 and Table 6, respectively. 

For both the MLCS and template fitting methods, the fit to the data determines the light curve 
parameters and their uncertainties. The "goodness" of the fits was within the expected statistical range with 
the exception of SN 1996J. This supernova is at a measured redshift of z=0.30, but some of the observations 
were obtained through a set of filters optimized for z=0.45. The uncertainty from this mismatch and the 
additional uncertainty from separate calibrations of the local standards' magnitudes in two sets of filters 
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may be the source of the poor result for this object. 

Four remaining SNe la presented here (SNe 1995ao, 1995ap, 1996T, and 1996R) are too sparsely 
sampled to provide meaningful light curves fit by cither of the light curve fitting methods. However, Riess et 
al. (1998a) describe a technique to measure the distance to sparsely observed SNe which lack well-sampled 
light curves. This "Snapshot" method measures the age and the luminosity/light-curve shape parameter 
from a SN la spectrum using techniques from Riess et al. (1997) and Nugent et al. (1995). An additional 
photometric epoch in 2 passbands (with host galaxy templates if needed) provides enough information 
to determine the extinction-free distance. For the four sparsely observed SNe la in our sample, we have 
measured the SN parameters with this method and list them in Table 7. This sample of sparsely observed, 
high-rcdshift SNe la is augmented by distances for SN 19971 (z = 0.17) and SN 1997ap (z = 0.83) from 
Riess et al. (1998a). 

For all SN la distance measurements, the dominant source of statistical uncertainty is the extinction 
measurement. The precision of our determination of the true extinction is improved using our prior 
understanding of its magnitude and direction (Riess, Press, & Kirshner 1996a, Appendix). 



4. Cosmological Implications of SNe la 
4.1. Cosmological Parameters 

Distance estimates from SN la light curves are derived from the luminosity distance, 

where £ and T are the SN's intrinsic luminosity and observed flux, respectively. In Fricdmann-Robertson- 
Walker cosmologies, the luminosity distance at a given rcdshift, z, is a function of the cosmological 
parameters. Limiting our consideration of these parameters to the Hubble constant, Ho, the mass density, 
Qm, and the vacuum energy density (i.e., the cosmological constant), Qa (but see Caldwell, Dave, & 
Stcinhardt 1998; Garnavich et al. 1998 for other energy densities), the luminosity distance is 



D L = cH^il + z) \n k \- 1/2 sinn{\VL k \ 1/2 f dz[(l + zf{\ + Sl M z) - z(2 + z)ft A ]~ 1/2 }, 

Jo 



(2) 



where Qfc = 1 — — ^A, and sinn is sinh for > and sin for Sl^ < (Carroll, Press, & Turner 1992). 
For Dl in units of Megaparsecs, the predicted distance modulus is 

fi p = 51ogD L + 25. (3) 

Using the data described in §2 and the fitting methods of §3 we have derived a set of distances, 
for SNe with 0.01 < z < 0.97. The available set of high-redshift SNe includes nine well-observed SNe la, 
six sparsely observed SNe la, and SN 1997ck (z — 0.97) whose light curve was well observed but lacks 
spectroscopic classification and color measurements. The Hubble diagrams for the nine well-observed SNe la 
plus SN 1997ck, with light curve distances calculated from the MLCS method and the template approach, 
are shown in Figures 4 and 5. The likelihood for the cosmological parameters can be determined from a \ 2 
statistic, where 

x 2 (^o, n M , o A ) = E M«-> H °^f?-»>J 2 ( 4 ) 
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and a v is the dispersion in galaxy redshift (in units of distance moduli) due to peculiar velocities. This 
term also includes the uncertainty in galaxy redshift. We have calculated this % 2 statistic for a wide range 
of the parameters H , Qm, and ^a- We do not consider the unphysical region of parameter space where 
CIm < 0; equation (2) describes the effect of massive particles on the luminosity distance. There is no reason 
to expect that the evaluation of equation (2) for VLm < has any correspondence to physical reality. We 
also neglect the region of (Qm,^a) parameter space which gives rise to so-called "bouncing" or rebounding 
Universes which do not monotonically expand from a "big bang" and for which equation (2) is not solvable 
(see Figure 6 and 7) (Carroll, Press, & Turner 1992). 

Due to the large redshifts of our distant sample and the abundance of objects in the nearby sample, our 
analysis is insensitive to a v within its likely range of 100 km s _1 < a v < 400 km s _1 (Marzke ct al. 1995; 
Lin et al. 1996). For our analysis we adopt a v — 200 km s _1 . For high-redshift SNe la whose redshifts were 
determined from the broad features in the SN spectrum (see Table 1), we add 2500 km s _1 in quadrature 
to a v . 

Separating the effects of matter density and vacuum energy density on the observed redshift-distance 
relation could in principle be accomplished with measurements of SNe la over a significant range of high 
redshifts (Goobar & Perlmutter 1995). Because the matter density decreases with time in an expanding 
Universe while the vacuum energy density remains constant, the relative influence of Qm to f^A on the 
redshift-distance relation is a function of redshift. The present data set has only a modest range of redshifts 
so we can only constrain specific cosmological models or regions of (12m, ^a) parameter space to useful 
precision. 

The x 2 statistic of equation (4) is well suited for determining the most likely values for the cosmological 
parameters Hq, Q,m, and Q A as well as the confidence intervals surrounding them. For constraining regions 
of parameter space not bounded by contours of uniform confidence (i.e., constant x 2 ), we need to define the 
probability density function (PDF) for the cosmological parameters. The PDF (p) of these parameters given 
our distance moduli is derived from the PDF of the distance moduli given our data from Bayes' theorem, 

p(H Q , n m , n A \fi ) = — , (5) 

where is our set of distance moduli (Lupton 1993). Since we have no prior constraints on the cosmological 
parameters (besides the excluded regions) or on the data, we take p(H , Q m , ^a) and p{no) to be constants. 
Thus we have for the allowed region of (H , Q m , ^a) 



^(iJc^m^AlMo) OC p(^ \H ,fl m ,fl A ). (6) 

We assume each distance modulus is independent (aside from systematic errors discussed in §5) and 
normally distributed, so the PDF for the set of distance moduli given the parameters is a product of 
Gaussians: 



p(Ho \H ,n m ,n A = I I =exp — - — — — . (7 



Ha 

Rewriting the product as a summation of the exponents and combining with equation (4) we have 
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p(fJ,o\H ,Q, m ,Q, A ) = \Y[ , 1 lexpf-^-Y (8) 

The product in front is a constant, so combining with equation (6) the PDF for the cosmological parameters 
yields the standard expression (Lupton 1993) 

p(HQ,£l m ,£l A \ii ) oc exp (^y) ' ^ 
The normalized PDF comes from dividing this relative PDF by its sum over all possible states, 



exp 



p(H , il m ,il A \p ) = — - — , (10) 

I-oo dH ° I-oo dfl A J CX P (-^-J dfl M 

neglecting the unphysical regions. The most likely values for the cosmological parameters and preferred 
regions of parameter space are located where equation (4) is minimized or alternately equation (10) is 
maximized. 

The Hubble constants as derived from the MLCS method, 65.2 ±1.3 km s" 1 Mpc -1 , and from the 
template fitting approach, 63.8 ±1.3 km s _1 Mpc^ 1 , are extremely robust and attest to the consistency 
of the methods. These determinations include only the statistical component of error resulting from the 
point-to-point variance of the measured Hubble flow and do not include any uncertainty in the absolute 
magnitude of SN la. From three photoclcctrically observed SNe la, SN 1972E, SN 1981B, and SN 1990N 
(Saha et al. 1994, 1997), the SN la absolute magnitude was calibrated from observations of Cepheids in 
the host galaxies. The calibration of the SN la magnitude from only three objects adds an additional 
5% uncertainty to the Hubble constant, independent of the uncertainty in the zeropoint of the distance 
scale. The uncertainty in the Cepheid distance scale adds an uncertainty of ~ 10% to the derived Hubble 
constant (Feast & Walker 1987; Kochanek 1997; Madore & Freedman 1998). A realistic determination of 
the Hubble constant from SNe la would give 65 ± 7 km s _1 Mpc -1 with the uncertainty dominated by the 
systematic uncertainties in the calibration of the SN la absolute magnitude. These determinations of the 
Hubble constant employ the Cepheid distance scale of Madore & Freedman (1991) which uses a distance 
modulus to the Large Magellanic Cloud (LMC) of 18.50 mag. Parallax measurements by the Hipparcos 
satellite indicate that the LMC distance could be greater and hence our inferred Hubble constant smaller 
by 5% to 10% (Reid 1997), though not all agree with the interpretation of these parallaxes (Madore & 
Freedman 1998). All subsequent indications in this paper for the cosmological parameters VLm and fl A are 
independent of the value for the Hubble constant or the calibration of the SN la absolute magnitude. 

Indications for Qm and Q A , independent from H , can be found by reducing our three-dimensional 
PDF to two dimensions. A joint confidence region for Q, M and Q A is derived from our three dimensional 
likelihood space 

/oo 
p(n M ,^A,H \fio)dH . (11) 
-oo 

The likelihood that the cosmological constant is greater than zero is given by summing the likelihood for 
this region of parameter space, 

p(n A > o\no) = / dn A / P (n M ,n A \no)dn M . (12) 
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This integral was evaluated numerically over a wide and finely spaced grid of cosmological parameters for 
which equation (11) is non-trivial. 

From the nine spectroscopic high-redshift SNe la with well-observed light and color curves, a 
non-negligible positive cosmological constant is strongly preferred at the 99.6% (2.9a) and >99.9% (3.9ct) 
confidence levels for the MLCS and template fitting methods, respectively (see Table 8). This region of 
parameter space is nearly identical to the one which results in an eternally expanding Universe. Boundless 
expansion occurs for a cosmological constant of 



Q A > 

(Carroll, Press, & Turner 1992), and its likelihood is 



< fl M < 1 

40 M {cos[i cos-^i^) + f ]} 3 n M > 1 



(13) 



pi pOO pOO pOO 

/ dn M p(fi M ,n A \^ Q )dn A + dn M p(n M ,n A \n )dn A . (w) 

JO JO Jl J4f2 M {cos[Icos- 1 (i 11 ^-) + ^ L ]} 3 

The preference for eternal expansion is numerically equivalent to the confidence levels cited for a 
non-negligible, positive cosmological constant. 

We can include external constraints on Q,m, ^a, or their sum to further refine our determination of the 
cosmological parameters. For a spatially flat Universe (i.e., Q.m + = ^tot = 1), we find £l A — 0.68 ± 0.10 
(fl M = 0.32 ± 0.10) and il A = 0.84 ± 0.09 (fl M = 0.16 ± 0.09) for MLCS and template fitting, respectively 
(see Table 8). The hypothesis that matter provides the closure density (i.e., Qm = 1) is ruled out at the la 
to 9ct level by either method. Again, f^A > and an eternally expanding Universe are strongly preferred, 
at this same confidence level. We emphasize that these constraints reflect statistical errors only; systematic 
uncertainties are confronted in §5. 

Other measurements based on the mass, light, x-ray emission, numbers, and motions of clusters of 
galaxies provide constraints on the mass density which have yielded typical values of Qm ~ 0.2 — 0.3 
(Carlberg et al. 1996; Bahcall, Fan, & Cen 1997; Lin et al. 1996; Strauss & Willick 1995). Using the 
constraint that 12m = 0.2 provides a significant indication for a cosmological constant: il A = 0.65 ± 0.22 
and f^A = 0.88 ±0.19 for the MLCS and template fitting methods, respectively (see Table 8). For Qm = 0.3 
we find il A = 0.80 ± 0.22 and il A = 0.96 ± 0.20 for the MLCS and template fitting methods, respectively. 

If we instead demand that fl A = 0, we are forced to relax the requirement that Qm > to locate 
a global minimum in our \ 2 statistic. Doing so yields an unphysical value of = —0.38 ± 0.22 and 
Qm = —0.52 ± 0.20 for the MLCS and template fitting approaches, respectively (see Table 8). This result 
emphasizes the need for a positive cosmological constant for a plausible fit. 

For the four sparsely observed SNe la (SN 1996R, SN 1996T, SN 1995ao, and SN 1995ap), we employed 
the snapshot distance method (Riess et al. 1998a) to determine the luminosity distances. Unfortunately, 
the low priority given to these objects resulted in observations not only limited in frequency but in 
signal-to-noise ratio as well. Consequently, these 4 distances are individually uncertain at the 0.4-0.6 mag 
level. We have compared these distances directly to a set of nine SNe la distances measured by the same 
snapshot method with 0.01 < z < 0.83 from Riess et al. (1998a) and reprinted here in Tables 7 and 9. This 
approach avoids the requirement that distances calculated from light curves and the snapshot method be 
on the same distance scale although this has been shown to be true (Riess et al. 1998a). 

The complete but sparse set of 13 snapshot distances now including six SNe la with z > 0.16 yields 
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conclusions which are less precise but fully consistent with the statistically independent results from the 
well-sampled SN la light curves (see Table 8). 

Having derived the two PDFs, p(^m,^a), f° r the ~ 40 SNc la light curves and the 13 incomplete 
( "snapshot" ) SNe la light curves independently, we can multiply the two PDFs to yield the PDF for all 
- 50 SNe la which includes 15 SNe with 0.16 < z < 0.62. Contours of constant PDF from the MLCS 
method and the template fitting method, each combined with the snapshot PDF, are shown in Figures 
6 and 7. These contours are closed by their intersection with the line SIm = and labeled by the total 
probability contained within. 

Including the snapshot distances modestly strengthens all of the previous conclusions about the 
detection of a non- negligible, positive cosmological constant (see Table 8). This set of 15 high-redshift SNe 
la favors J7a > and an eternally expanding Universe at 99.7% (3.0a) and >99.9% (4. Oct) confidence for 
the MLCS and template fitting methods, respectively. This complete set of spectroscopic SNe la represents 
the full strength of the high-redshift sample and provides the most reliable results. 

A remarkably high-redshift supernova (z = 0.97), SN 1997ck, was excluded from all these analyses due 
to its uncertain extinction and the absence of a spectroscopic identification. Nevertheless, if we assume a 
negligible extinction of Ab = 0.0 ± 0.1 for SN 1997ck as observed for the rest of our high-redshift sample 
and further assume it is of type la, as its well-observed B rest-frame light curve suggests (see Figure 
3), we could include this object in our previous analysis (see Tabic 8). As seen in Figures 6 and 7, SN 
1997ck constrains specific values of Q.m and f^A by effectively closing our confidence contours because of 
the increased redshift range of this augmented sample. The values implied using SN 1997ck and the rest of 
the spectroscopic SNe la, under the previous assumptions, are f^M = 0.241^24, £1a = 0.72^48 from the 
MLCS method and Qm — O^Oi^g, CI a = 1-56+q;7q from the template fitting method. The preference for a 
non-negligible, positive cosmological constant remains strong (see Table 8). 

As seen in Table 8, the values of the xt f° r the cosmological fits are reassuringly close to unity. This 
statement is more meaningful for the MLCS distances which are accompanied by statistically reliable 
estimates of the distance uncertainty (Riess, Press, & Kirshner 1996a). The values for xt indicate a good 
agreement between the expected distance uncertainties and the observed distance dispersions around the 
best fit model. They leave little room for sources of additional variance, as might be introduced by a 
significant difference between the properties of SN la at high and low redshift. 

4.2. Deceleration Parameter 

An alternate approach to exploring the expansion history of the Universe is to measure the current 
(z = 0) deceleration parameter, q = —a(t )a(to)/a 2 (to), where a is the cosmic scale factor. Because 
the deceleration is defined at the current epoch and the supernovae in our sample cover a wide range in 
redshift, we can only determine the value of qo within the context of a model for its origin. Nevertheless, 
for moderate values of deceleration (or acceleration) the determination of qo from our SNe, all of which are 
at z < 1, provides a valuable description of the current deceleration parameter valid for most equations of 
state of the Universe. 

We have derived estimates of q within a two-component model where q = ^f- — £1a- This definition 
assumes that the only sources of the current deceleration are mass density and the cosmological constant. 
A more complete definition for qo would include all possible forms of energy density (see Caldwell, Dave, & 
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Stcinhardt 1998) but is beyond the scope of this paper. From our working definition of qo, negative values 
for the current deceleration (i.e., accelerations) are generated only by a positive cosmological constant and 
not from unphysical, negative mass density. 

Current acceleration of the expansion occurs for a cosmological constant of 

> (15) 

/•oo />oo 

P(q < 0|/i o ) = / dQ M / p(f7 M ,r2 A |Ato)e^A (16) 
Jo J^- 



and its likelihood is 



considering only Qm > 0. Figures 6 and 7 show the boundary between current acceleration and deceleration 
as well as lines of constant go- Fot the complete set of supernova distances (excluding SN 1997ck), current 
acceleration is strongly preferred at the 99.5% (2.8a) confidence level for the MLCS method and >99.9% 
level (3.9a) for the template fitting approach. The most likely value for q is given by the peak of the 
distribution 

poo poo 

p(Qo\^o) = / dn M / p(fiM,fU|go,Mo)<2fU- (17) 

Jo J -oo 

This expression determines the likelihood of a given q from the sum of the likelihoods of the combinations 
of fijvf and which produce that value of q . Values for q and their uncertainties for the different 
methods and sample cuts are summarized in Table 8. With the current sample we find a robust indication 
for the sign of qo and a more uncertain estimate for its value, qo = — 1.0 ± 0.4. Because lines of constant 
qo are skewed with respect to the major axis of our uncertainty contours, more SNe la at redshifts greater 
than z = 0.5 will be needed to yield a more robust indication for the value of qo. 



4.3. Dynamical Age of the Universe 

The dynamical age of the Universe can be calculated from the cosmological parameters. In an empty 
Universe with no cosmological constant, the dynamical age is simply the inverse of the Hubble constant; 
there is no deceleration. SNe la have been used to map the nearby Hubble flow resulting in a precise 
determination of the Hubble constant (Hamuy et al. 1995, 1996a; Riess, Press, & Kirshner 1995, 1996a). 
For a more complex cosmology, integrating the velocity of the expansion from the current epoch (z = 0) to 
the beginning (z = oo) yields an expression for the dynamical age 

f-OO 

t (H ,n M ,^A) = H^ 1 (l + z)- 1 [(l + zf(l + n M z)-z(2 + z)n A ]- 1 / 2 dz (18) 
Jo 

(Carroll, Press, & Turner 1992). Combining a PDF for the cosmological parameters, p(Ho,^m,^a), with 
the above expression we can derive the PDF for the age of the Universe: 

/oo poo P oo 

dH / dCt M / p(H ,Q M ,^A\to,fJi})d^A- (19) 
-oo JO J —oo 

Equation (19) expresses the likelihood for a given age, t 0} as the sum of the likelihoods of all combinations 
of Ho, CIm, and Q\ which result in the given age. The peak of this function provides our maximum 
likelihood estimate for the dynamical age, to- Without SN 1997ck, the peak is at 13.6^0.8 Gyr from the 
MLCS PDF. For the template fitting approach the peak occurs at 14.8lJ;g Gyr. A naive combination of the 
two distributions yields an estimate of 14.2+Q g Gyr adopting either method's uncertainty (see Figure 8). 
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Again, these errors include only the statistical uncertainties of the measurement. Including the systematic 
uncertainty of the Cepheid distance scale, which may be as much as 10%, a reasonable estimate of the 
dynamical age would be 14.2±1.5 Gyr. 

An illuminating way to characterize the dynamical age independent of the Hubble constant is to 
measure the product H t - For the MLCS method, the template fitting method, and the combination of the 
two, we find H to to be 0.90, 0.96, and 0.93, respectively. These values imply a substantially older Universe 
for a given value of Hq in better accordance with globular cluster ages than the canonical value of H^t^=2/'i 
for Q.m = 1 and Cl\ = 0. Our determination of the dynamical age of the Universe is consistent with the 
rather wide range of values of the ages using stellar theory or radioactive dating. Oswalt et al. (1996) have 
shown that the Galactic disk has a lower age limit of 9.5 Gyr measured from the cooling sequence of the 
white dwarfs. The radioactive dating of stars via the thorium and europium abundances gives a value of 
15.2 ± 3.7 Gyr (Cowan et al. 1997). We can expect these ages to become more precise as more objects are 
observed. 

Perhaps the most widely quoted ages of the Universe come from the age estimates of globular cluster 
stars. These are dependent on the distance scale used and the stellar models employed. Vandenberg, 
Stetson, & Bolte (1996) note that these two effects generally work in the opposite direction: for instance, 
if one increases the distance to the Large Magellanic Cloud (LMC), the dynamical age of the Universe 
increases while the age based on the cluster ages decreases (the main-sequence turnoff is brighter implying 
a younger population). This means that there is only a limited range in cosmological and stellar models 
that can bring the two ages into concordance. 

Prior to Hipparcos, typical age estimates based on the subdwarf distance scale were greater than 15 Gyr 
for cluster ages. Bolte & Hogan (1995) find 15.8 ±2.1 Gyr for the ages of the oldest clusters, while Chaboyer, 
Dcmarque, & Sarajedini (1996) find a typical age of 18 Gyr for the oldest clusters. Chaboyer (1995) also 
estimates the full range of viable ages to be 11-21 Gyr with the dominant error due to uncertainties in the 
theory of convection. An independent distance scale based on parallaxes of white dwarfs provides an age 
estimate for the globular cluster M4 of 14.5 — 15.5 Gyr (Renzini et al. 1996). 

However, the Hipparcos parallaxes suggest an increased distance to the LMC and the globular clusters 
(Feast & Catchpole 1997; Reid 1997; but see Madore & Freedman 1998). With this new distance scale, the 
ages of the clusters have decreased to about 11.5 Gyr with an uncertainty of 2 Gyr (Gratton et al. 1997; 
Chaboyer et al. 1998). Given the large range in ages from the theoretical models of cluster turnoffs and 
the inconsistency of the subdwarf and white dwarf distance scales applied to the ages of globular clusters, 
a robust estimate for the ages of the globular clusters remains elusive. Even with these uncertainties, the 
dynamical age of the Universe derived here is consistent with the ages based on stellar theory or radioactive 
dating. Evidently, there is no longer a problem that the age of the oldest stars is greater than the dynamical 
age of the Universe. 

Despite our inability to place strong constraints on the values for Qm and f^A independently, our 
experiment is sensitive to the difference of these parameters. Because the dynamical age also varies 
approximately as the difference in Qm and Q\, our leverage on the determination of the dynamical age is 
substantial. This point can be illustrated with a display of lines of constant dynamical age as a function of 
Q.M and comparing Figure 9 to Figures 6 and 7, we see that the semi-major axes of our error ellipses 
are roughly parallel to the lines of constant dynamical age. Figure 9 also indicates why the most likely value 
for the dynamical age differs from the dynamical age derived for the most likely values of Hq, I^m, and 
JIa- For a fixed value of the Hubble constant, younger dynamical ages span a larger region of the (Qm^a) 
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parameter space than older ages. This shifts the most likely value for to towards a younger age and results 
in a "tail" in the distribution, p(t \fi ), extending towards older ages. 

5. Discussion 

The results of §4 suggest an eternally expanding Universe which is accelerated by energy in the 
vacuum. Although these data do not provide independent constraints on Q,m and J7a to high precision 
without ancillary assumptions or inclusion of a supernova with uncertain classification, specific cosmological 
scenarios can still be tested without these requirements. 

High-redshift SNe la are observed to be dimmer than expected in an empty Universe (i.e., Q,m = 0) 
with no cosmological constant. A cosmological explanation for this observation is that a positive vacuum 
energy density accelerates the expansion. Mass density in the Universe exacerbates this problem, requiring 
even more vacuum energy. For a Universe with VIm = 0.2, the MLCS and template fitting distances to the 
well-observed SNe are 0.25 and 0.28 mag farther on average than the prediction from J7a = 0. The average 
MLCS and template fitting distances are still 0.18 and 0.23 mag farther than required for a 68.3% (la) 
consistency for a Universe with Qm = 0.2 and without a cosmological constant. 

Depending on the method used to measure all the spectroscopically confirmed SN la distances, we find 
Qa to be inconsistent with zero at the 99.7% (3. Oct) to >99.9% (4. Oct) confidence level. Current acceleration 
of the expansion is preferred at the 99.5% (2.8a) to >99.9% (3.9ct) confidence level. The ultimate fate of the 
Universe is sealed by a positive cosmological constant. Without a restoring force provided by a surprisingly 
large mass density (i.e., Qm > 1) the Universe will continue to expand forever. 

How reliable is this conclusion? Although the statistical inference is strong, here we explore systematic 
uncertainties in our results with special attention to those that can lead to overestimates of the SNe la 
distances. 



5.1. Evolution 

The local sample of SNe la displays a weak correlation between light curve shape (or luminosity) and 
host galaxy type. The sense of the correlation is that the most luminous SNe la with the broadest light 
curves only occur in late-type galaxies. Both early-type and late-type galaxies provide hosts for dimmer SNe 
la with narrower light curves (Hamuy et al. 1996c). The mean luminosity difference for SNe la in late- type 
and early-type galaxies is ~ 0.3 mag (Hamuy et al. 1996c). In addition, the SN la rate per unit luminosity 
is almost twice as high in late-type galaxies as in early-type galaxies at the present epoch (Cappellaro et 
al. 1997). This suggests that a population of progenitors may exist in late-type galaxies which is younger 
and gives rise to brighter SNe la (with broader light curves) than those contained in early-type galaxies or 
within pockets of an older stellar population in the late-type galaxies. Such observations could indicate an 
evolution of SNe la with progenitor age. 

Hoflich, Thielemann, & Wheeler (1998) calculate differences in the light curve shape, luminosity, and 
spectral characteristics of SNe la as a function of the initial composition and metallicity of the white 
dwarf progenitor. As we observe more distant samples, we expect the progenitors of SN la to come 
from a younger and more metal-poor population of stars. Hoflich, Thielemann, & Wheeler (1998) have 
shown that a reduction in progenitor metallicity by a factor of 3 has little effect on the SN la bolometric 
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luminosity at maximum. For their models, such a change in metallicity can alter the peak luminosity by 
small amounts (~ 0.05 mag) in rest-frame B and V, accompanied by detectable spectral signatures. These 
spectral indicators of evolution are expected to be most discernible in the rest-frame U passband where 
line blanketing is prevalent. Future detailed spectral analyses at these short wavelengths might provide a 
constraint on a variation in progenitor metallicity. 

The effect of a decrease in SN la progenitor age at high rcdshift is predicted to be more significant 
than metallicity (Hoflich, Thielemann, & Wheeler 1998). Younger white dwarfs are expected to evolve from 
more massive stars with a lower ratio of C/O in their cores. The lower C/O ratio of the white dwarf reduces 
the amount of 56 Ni synthesized in the explosion, but an anticipated slower rise to maximum conserves more 
energy for an increased maximum brightness. By reducing the C/O ratio from 1/1 to 2/3, the B — V color 
at maximum is expected to become redder by 0.02 mag and the post-maximum decline would become 
steeper. This prediction of a brighter SN la exhibiting a faster post-maximum decline is opposite to what 
is seen in the nearby sample (Phillips 1993; Hamuy et al. 1995; Hamuy et al. 1996a, b,c,d; Riess, Press, & 
Kirshner 1996a; Appendix) and will be readily testable for an enlarged high redshift sample. Specifically, 
a larger sample of distant SNe la (currently being compiled) would allow us to determine the light curve 
shape relations at high-redshift and test whether these evolve with look-back time. Presently, our sample is 
to small to make such a test meaningful. 

We expect that the relation between light curve shape and luminosity that applies to the range of 
stellar populations and progenitor ages encountered in the late-type and early-type hosts in our nearby 
sample should also be applicable to the range we encounter in our distant sample. In fact, the range of 
age for SN la progenitors in the nearby sample is likely to be larger than the change in mean progenitor 
age over the 4 to 6 Gyr look-back time to the high-redshift sample. Thus, to first order at least, our local 
sample should correct our distances for progenitor or age effects. 

We can place empirical constraints on the effect that a change in the progenitor age would have on 
our SN la distances by comparing subsamples of low redshift SNe la believed to arise from old and young 
progenitors. In the nearby sample, the mean difference between the distances for the early-type (8 SNe 
la) and late-type hosts (19 SNe la), at a given redshift, is 0.04 ± 0.07 mag from the MLCS method. This 
difference is consistent with zero. Even if the SN la progenitors evolved from one population at low redshift 
to the other at high redshift, we still would not explain the surplus in mean distance of 0.25 mag over 
the f^A = prediction. For the template fitting approach, the mean difference in distance for SNe la in 
early-type and late-type hosts is 0.05 ± 0.07 mag. Again, evolution provides an inadequate explanation for 
the 0.28 mag difference in the template fitting SNe la distances and the fl\ = prediction. 

However, the low-redshift sample is dominated by late-type hosts and these may contain a number of 
older progenitors. It is therefore difficult to assess the precise effect of a decrease in progenitor age at high 
redshift from the consistency of distances to early- type and late-type hosts (see Schmidt et al. 1998). If, 
however, we believed that young progenitors give rise to brighter SNe la with broader light curves (Hamuy 
et al. 1996c) as discussed above, we could more directly determine the effect on distance determinations 
of drawing our high-redshift sample from an increasingly youthful population of progenitors. The mean 
difference in the Hubble line defined by the full nearby sample and the subsample of SNe la with broader 
than typical light curves (A < 0) is 0.02 ± 0.07 for the MLCS method. For the template fitting method, 
the difference between the full sample and those with broader light curves (Ami 5 (B) < 1.1) is 0.07 ± 
0.07. Again, we find no indication of a systematic change in our distance estimates with a property that 
may correspond to a decrease in progenitor age. Another valuable test would be to compare low-redshift 
distances to starburst and irregular type galaxies which presumably are hosts to progenitors which are 
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young and metal-poor. Such a nearby sample may yield the closest approximation to the SNe la observed 
at high redshift. Future work will be needed to gather this informative sample which would be composed 
of objects such as SN 1972E in NGC 5253 (which we does fit the luminosity light curve shape relations; 
Hamuy et al 1996b). 

Another check on evolutionary effects is to test whether the distribution of light curve decline rates 
is similar between the nearby sample of supernovae and the high-redshift sample. Figure 10 shows the 
observed distribution of the MLCS light-curve shape parameters, A, and the template fitting parameters, 
Ami5(_B), with redshift. A Kolmogorov-Smirnov test shows no significant difference in the distributions of 
the low and high-redshift samples, but the sample is too small to be statistically significant. The actual 
difference in mean luminosity between the low-redshift and high-redshift samples implied by the light curve 
shapes is 0.02 mag by either method. We conclude that there is no obvious difference between the shapes 
of SNe la light curves at z « and at z « 0.5. 

It is reassuring that initial comparisons of high-redshift SN la spectra appear remarkably similar to 
those observed at low-redshift. This can be seen in the high signal-to-noise ratio spectra of SN 1995ao 
(z = 0.30) and SN 1995ap (z — 0.23) in Figure 1. Another demonstration of this similarity at even higher 
redshift is shown in Figure 11 for SN 1998ai (z = 0.49; IAUC 6861) whose light curve was not used in this 
work. The spectrum of SN 1998ai was obtained at the Keck telescope with a 5 x 1800 s exposure using 
LRIS and was reduced as described in §2.2 (Filippenko et al. 1998). The spectral characteristics of this 
SN la appear to be indistinguishable from the range of characteristics at low redshift to good precision. In 
additon, a time sequence of spectra of SN la 1997ex (z=0.36; Nugent et al. 1998a) compared with those of 
local SNe la reveals no significant spectral differences (Filippenko et al. 1998). 

We expect that our local calibration will work well at eliminating any pernicious drift in the supernova 
distances between the local and distant samples. Until we know more about the stellar ancestors of SNe 
la, we need to be vigilant for changes in the properties of the supernovae at significant look-back times. 
Our distance measurements could be particularly sensitive to changes in the colors of SNe la for a given 
light curve shape. Although our current observations reveal no indication of evolution of SNe la at z « 0.5, 
evolution remains a serious concern which can only be eased and perhaps understood by future studies. 

5.2. Extinction 

Our SNe la distances have the important advantage of including corrections for interstellar extinction 
occurring in the host galaxy and the Milky Way. The uncertainty in the extinctions is a significant 
component of error in our distance uncertainties. Extinction corrections based on the relation between 
SN la colors and luminosity improve distance precision for a sample of SNe la that includes objects with 
substantial extinction (Riess, Press, & Kirshner 1996a). Yet, in practice, we have found negligible extinction 
to the high-redshift SNe la. The mean B — V color at maximum is —0.13 ± 0.05 from the MLCS method 
and —0.07 ± 0.05 from the template fitting approach, consistent with an unreddened B — V color of —0.10 
to —0.05 expected for slowly declining light curves as observed in the high-redshift sample (Riess, Press, & 
Kirshner 1996a; Appendix). 

Further, the consistency of the measured Hubble flow from SNe la with late-type and early-type 
hosts (§5.1) shows that the extinction corrections applied to dusty SNe la at low redshift do not alter the 
expansion rate from its value measured from SNe la in low dust environments. The conclusions reached 
in §4 would not alter if low and high-redshift SNe with significant extinction were discarded rather than 



- 21 - 



included after a correction for extinction. 

The results of §4 do not depend on the value of the ratios between color excess and selective absorption 
used to determine the extinctions of the high-redshift sample because the mean observed reddening is 
negligible. Some modest departures from the Galactic reddening ratios have been observed in the Small 
and Large Magellanic Clouds, M31, and the Galaxy, and they have been linked to metallicity variations 
(Walterbos 1986; Hodge & Kennicutt 1982; Bouchet et al 1985; Savage & Mathis 1979). Although our 
current understanding of the reddening ratios of interstellar dust at high redshift is limited, the lack of 
any significant color excess observed in the high-redshift sample indicates that the type of interstellar dust 
which reddens optical light is not obscuring our view of these objects. 

Riess, Press, & Kirshner (1996b) found indications that the Galactic ratios between selective absorption 
and color excess are similar for host galaxies in the nearby (z < 0.1) Hubble flow. Yet, what if these ratios 
changed with look-back time? Could an evolution in dust grain size descending from ancestral interstellar 
"pebbles" at higher redshifts cause us to underestimate the extinction? Large dust grains would not imprint 
the reddening signature of typical interstellar extinction upon which our corrections rely. However, viewing 
our SNe through such grey interstellar grains would also induce a dispersion in the derived distances. To 
estimate the size of the dispersion, we assume that the grey extinction is distributed in galaxies in the same 
way as typical interstellar extinction. 

Hatano, Branch, & Deaton (1997) have calculated the expected distribution of SN la extinction along 
random lines of sight in the host galaxies. A grey extinction distribution similar to theirs could yield 
differing amounts of mean grey extinction depending on the likelihood assigned to observing an extinction 
of Ab=0.0 mag. In the following calculations we vary only the likelihood of Ab— 0.0 mag to derive new 
extinction distributions with varying means. These different distributions also have differing dispersions of 
extinction. A mean grey extinction of 0.25 mag would be required to explain the measured MLCS distances 
without a cosmological constant. Yet the dispersion of individual extinctions for a distribution with a 
mean of 0.25 mag would be er^^O.40 mag, significantly larger than the 0.21 mag dispersion observed in 
the high-redshift MLCS distances. Grey extinction is an even less likely culprit with the template fitting 
approach; a distribution with a mean grey extinction of 0.28 mag, needed to replace a cosmological constant, 
would yield a dispersion of 0.42 mag, significantly higher than the distance dispersion of 0.17 mag observed 
in the high-redshift template fitting distances. 

Furthermore, most of the observed scatter is already consistent with the estimated statistical errors as 
evidenced by the xl (Table 8), leaving little to be caused by grey extinction. Nevertheless, if we assumed 
that all of the observed scatter were due to grey extinction, the mean shift in the SNe la distances would 
only be 0.05 mag. With the observations presented here, we cannot rule out this modest amount of grey 
interstellar extinction. 

This argument applies not only to exotic grey extinction but to any interstellar extinction not accounted 
for which obscures SNe la. Any spotty interstellar extinction which varies with line-of-sight in a way similar 
to the Hatano, Branch, & Deaton (1997) model of galaxies will add dispersion to the SN la distances. The 
low dispersion measured for the high-redshift sample places a strong limit on any mean spotty interstellar 
extinction. 

Grey intergalactic extinction could dim the SNe without either telltale reddening or dispersion, if all 
lines of sight to a given redshift had a similar column density of absorbing material. The component of the 
intergalactic medium with such uniform coverage corresponds to the gas clouds producing Lyman-a forest 
absorption at low redshifts. These clouds have individual H I column densities less than about 10 15 cm~ 2 
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(Bahcall et al. 1996). However, these clouds display low metallicities, typically less than 10% of solar. Grey 
extinction would require larger dust grains which would need a larger mass in heavy elements than typical 
interstellar grain size distributions to achieve a given extinction. Furthermore, these clouds reside in hard 
radiation environments hostile to the survival of dust grains. Finally, the existence of grey intergalactic 
extinction would only augment the already surprising excess of galaxies in high-redshift galaxy surveys 
(Huang et al. 1997). 

We conclude that grey extinction does not seem to provide an observationally or physically plausible 
explanation for the observed faintness of high-redshift SNe la. 

5.3. Selection Bias 

Sample selection has the potential to distort the comparison of nearby and distant supernovae. Most 
of our nearby (z < 0.1) sample of SNe la was gathered from the Calan/Tololo survey (Hamuy et al. 1993a) 
which employed the blinking of photographic plates obtained at different epochs with Schmidt telescopes 
and from less well-defined searches (Riess et al. 1998b). Our distant (z > 0.16) sample was obtained by 
subtracting digital CCD images at different epochs with the same instrument setup. 

If they were limited by the flux of the detected events, both nearby and distant SN la searches would 
preferentially select intrinsically luminous objects because of the larger volume of space in which these 
objects can be detected. This well-understood selection effect could be further complicated by the properties 
of SNe la; more luminous supernovae have broader light curves (Phillips 1993; Hamuy et al. 1995, 1996c; 
Riess, Press, & Kirshner 1995, 1996a). The brighter supernovae remain above a detection limit longer 
than their fainter siblings, yet also can fail to rise above the detection limit in the time interval between 
successive search epochs. The complex process by which SNe la are selected in low and high-redshift 
searches can be best understood with simulations (Hamuy & Pinto 1998). Although selection effects could 
alter the ratio of intrinsically dim to bright SNe la in our samples relative to the true population, our use 
of the light curve shape to determine the supernova's luminosity should correct most of this selection bias 
on our distance estimates. However, even after our light-curve shape correction, SNe la still have a small 
dispersion as distance indicators (cr « 0.15 mag), and any search program would still preferentially select 
objects which are brighter than average for a particular light curve shape and possibly select objects whose 
light curve shapes aid detection. 

To investigate the consequence of sample selection effects, we used a Monte Carlo simulation to 
understand how SNe la in our nearby and distant samples were chosen. For the purpose of this simulation 
we first assumed that the SN la rate is constant with look-back time. We assembled a population of SNe la 
with luminosities described by a Gaussian random variable um b = 0.4 mag and light-curve shapes which 
correspond to these luminosities as described by the MLCS vectors (see the Appendix). A Gaussian random 
uncertainty of a = 0.15 mag is assumed in the determination of absolute magnitude from the shape of a 
supernova's light curve. The time interval between successive search epochs, the search epoch's limiting 
magnitudes, and the apparent light-curve shapes were used to determine which SNe la were "discovered" 
and included in the simulation sample. A separate simulation was used to select nearby objects, with the 
appropriate time interval between epochs and estimates of limiting magnitudes. The results are extremely 
encouraging, with recovered values exceeding the simulated value of VLm or £1a by only 0.02 for these two 
parameters considered separately. Smoothly increasing the SN la rate by a factor of 10 by z = 1 doubles 
this bias to 0.04 for either parameter. There are two reasons we find such a small selection bias in the 
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recovered cosmological parameters. First, the small dispersion of our distance indicator results in only a 
modest selection bias. Second, both nearby and distant samples include an excess of brighter than average 
SNe, so the difference in their individual selection biases remains small. 

As discussed by Schmidt et al. (1998), obtaining accurate limiting magnitudes is complex for the 
CCD-based searches, and essentially impossible for the photographic searches. Limiting magnitudes vary 
from frame to frame, night to night, and film to film, so it is difficult to use the actual detection limits in our 
simulation. Nevertheless, we have run simulations varying the limiting magnitude, and this does not change 
the results significantly. We have also tried increasing the dispersion in the SN la light curve shape vs. 
absolute magnitude correlation at wavelengths shorter than 5000 A. Even doubling the distance dispersion 
of SNe la (as may be the case for rest-frame U) does not significantly change the simulation results. 

Although these simulations bode well for using SNe la to measure cosmological parameters, there 
are other differences between the way nearby and distant supernova samples are selected which are more 
difficult to model and are not included in our present simulations. Von Hippel, Bothun, & Schommer 
(1997) have shown that the selection function of the nearby searches is not consistent with that of a 
strict magnitude-limited search. It is unclear whether a photographic search selects SNe la with different 
parameters or environments than a CCD search or how this could affect a comparison of samples. Future 
work on quantifying the selection criteria of the samples is needed. A CCD search for SNe la in Abell 
clusters by Reiss et al. (1998) will soon provide a nearby SN la sample with better understood selection 
criteria. Although indications from the distributions of SN la parameters suggest that both our searches 
have sampled the same underlying population (see Figure 10), we must continue to be wary of subtle 
selection effects which might bias the comparison of SNe la near and far. 



5.4. Effect of a Local Void 

It has been noted by Zehavi et al. (1998) that the SNe la out to 7000 km s _1 exhibit an expansion 
rate which is 6% greater than that measured for the more distant objects. The significance of this peculiar 
monopolc is at the 2a to 3a confidence level; it is not inconsistent with the upper limit of ~ 10% for the 
difference between the local and global values of H found by Kim et al. (1997). The implication is that 
the volume out to this distance is underdense relative to the global mean density. This effect appears as an 
excess redshift for a given distance modulus (within 7000 km s _1 ) and can be seen with both the MLCS 
method and the template fitting method in Figures 4 and 5 . 

If true, what effect would this result have on our conclusions? In principle, a local void would 
increase the expansion rate measured for our low-redshift sample relative to the true, global expansion 
rate. Mistaking this inflated rate for the global value would give the false impression of an increase in the 
low-redshift expansion rate relative to the high-rcdshift expansion rate. This outcome could be incorrectly 
attributed to the influence of a positive cosmological constant. In practice, only a small fraction of our 
nearby sample is within this local void, reducing its effect on the determination of the low-redshift expansion 
rate. 

As a test of the effect of a local void on our constraints for the cosmological parameters, we reanalyzed 
the data discarding the seven SNe la within 7000 km s _1 (108 Mpc for H — 65). The result was a 
reduction in the confidence that fi A > from 99.7% (3. Oct) to 98.3% (2.4<r) for the MLCS method and from 
>99.9% (4.0a) to 99.8% (3.1cr) for the template fitting approach. The tests for both methods excluded the 
unclassified SN 1997ck and included the snapshot sample, the latter without two SNe la within 7000 km 
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s 1 . As expected, the influence of a possible local void on our cosmological conclusions is relatively small. 

5.5. Weak Gravitational Lensing 

The magnification and demagnification of light by large-scale structure can alter the observed 
magnitudes of high-rcdshift supernovae (Kantowski, Vaughan, & Branch 1995). The effect of weak 
gravitational lensing on our analysis has been quantified by Wambsganss et al. (1997) and summarized by 
Schmidt et al. (1998). SN la light will, on average, be demagnified by 0.5% at z = 0.5 and 1% at z = 1 in 
a Universe with a non- negligible cosmological constant. Although the sign of the effect is the same as the 
influence of a cosmological constant, the size of the effect is negligible. 

Holz & Wald (1997) have calculated the weak lensing effects on supernova light from ordinary matter 
which is not smoothly distributed in galaxies but rather clumped into stars (i.e., dark matter contained in 
MACHOS). With this scenario, microlensing by compact masses becomes a more important effect further 
decreasing the observed supernova luminosities at z = 0.5 by 0.02 mag for fiM=0.2 (Holz 1998). Even if 
most ordinary matter were contained in compact objects, this effect would not be large enough to reconcile 
the SNe la distances with the influence of ordinary matter alone. 

5.6. Light Curve Fitting Method 

As described in §3.2, two different light curve fitting methods, MLCS (Riess, Press, & Kirshncr 1996a; 
Appendix) and a template fitting approach (Hamuy et al. 1995, 1996d), were employed to determine the 
distances to the nearby and high-redshift samples. Both methods use relations between light curve shape 
and luminosity as determined from SNe la in the nearby Hubble flow. Both methods employ an extinction 
correction from the measured color excess using relations between intrinsic color and light curve shape. 
In addition, both the MLCS and template fitting methods yield highly consistent measurements for the 
Hubble constant of i?o=65.2 ±1.3 and Hq—63.8 ±1.3, respectively not including any uncertainty in the 
determination of the SN la absolute magnitude which is the dominant uncertainty. It is also worth noting 
that both methods yield SN la distance dispersions of ~ 0.15 mag when complete light curves in B,V,R, 
and I are employed. For the purpose of comparing the same data at high and low rcdshifts, the use of SN 
la observations at low redshift were restricted to only B and V within 40 days of maximum light. 

Although the conclusions reached by the two methods when applied to the high-redshift SNe arc highly 
consistent, some differences are worth noting. There are small differences in the distance predictions at high 
redshift. For the distant sample, the template fitting distances exhibit a scatter of 0.17 mag around the 
best fit model as compared to 0.21 mag for the MLCS method. In addition, the template fitting distances 
to the high-redshift SNe la are on (weighted) average 0.03 mag farther than the MLCS distances relative 
to the low-redshift sample. These differences together result in slightly different confidence intervals for the 
two methods (see Figures 6, 7, and 8 and Table 8). For the set of 10 well-observed SNe la, a sample with 
scatter 0.17 mag or less is drawn from a population of scatter 0.21 mag 25% of the time. The chance that 
10 objects could be drawn from this same population with a mean difference of 0.03 mag is 66%. Future 
samples of SNe la will reveal if the observed differences arc explained by chance. Until then, we must 
consider the difference between the cosmological constraints reached from the two fitting methods to be a 
systematic uncertainty. Yet, for the data considered here, both distance fitting methods unanimously favor 
the existence of a non-negligible, positive cosmological constant and an accelerating Universe. 
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5.7. Sample Contamination 

The mean brightness of SNe la is typically 4 to 40 times greater than that of any other type of 
supernova, favoring their detection in the volume of space searched at high redshift. Yet in the course of 
our high-rcdshift supernova search (and that of the Supernova Cosmology Project; Perlmutter et al. 1995) 
a small minority of other supernova types have been found and we must be careful not to include such 
objects in our SN la sample. The classification of a supernova is determined from the presence or absence 
of specific features in the spectrum (Wheeler & Harkness 1990; Branch, Fisher, & Nugent 1993; Filippcnko 
1997). The spectra of Type la supernovae show broad Si II absorption near 6150 A Ca II (H&K) absorption 
near 3800 A a S II absorption doublet near 5300 A and 5500 A and numerous other absorption features 
with ionized Fe a major contributor (Filippcnko 1997). For supernovae at high redshift, some of these 
characteristic features shift out of the observer's frequency range as other, shorter wavelength features 
become visible. Classification is further complicated by low signal-to-noise ratio in the spectra of distant 
objects. The spectra of SNe la evolve with time along a remarkably reliable sequence (Riess et al. 1997). 
Final spectral classification is optimized by comparing the observed spectrum to well-observed spectra of 
SNe la at the same age as determined from the light curves. 

For most of the spectra in Figure 1, the identification as a SN la is unambiguous. However, in three 
of the lowest signal-to-noise ratio cases - 1996E, 1996H, and 19961 - the wavelengths near Si II absorption 
(rest-frame 6150 A ) were poorly observed and their classification warrants closer scrutiny. These spectra 
are inconsistent with Type II spectra which show H/3 (4861 A ) in emission and absorption and lack Fe 
II features shortly after maximum. These spectra are also inconsistent with Type lb spectra which would 
display He I A5876 absorption at a rest wavelength of ~ 5700 A. 

The most likely supernova type to be misconstrued as a Type la is a Type Ic, as this type comes closest 
to matching the SN la spectral characteristics. Although SN Ic spectra lack Si II and S II absorption, 
the maximum-light spectra at blue wavelengths can resemble those of SNe la ~ 2 weeks past maximum 
when both are dominated by absorption lines of Fe II with P Cygni profiles. Type Ic events are rare and 
one luminous enough to be found in our search would be rare indeed, but not without precedent. An 
example of such an object is SN 1992ar (Clocchiatti et al. 1998), which was discovered in the course of 
the Calan/Tololo SN survey and which reached an absolute magnitude, uncorrected for host galaxy dust 
extinction, of M v = -19.3 (H = 65 km s" 1 Mpc" 1 ). For both SN 1996H and SN 19961, the spectral 
match with a Type la at rest wavelengths less than 4500 A is superior to the fit to a Type Ic spectrum 
(see Figure 1). In both cases the spectra rise from deep troughs at the 3800 A Ca II break (rest-frame) to 
strong peaks at 3900 to 4100 A (rest-frame) as observed in SNe la. Type Ic spectra, by comparison, tend 
to exhibit a much weaker transition from trough to peak redward of the Ca II break (see Figure 12). 

For SN 1996E, the spectral coverage does not extend blueward of a rest wavelength of 4225 A rendering 
this diagnostic unusable. The absence of pre-maximum observations of SN 1996E makes it difficult to 
determine the age of the spectrum and that of the appropriate comparison spectra. As shown in Figure 12, 
the spectroscopic and photometric data for SN 1996E are consistent with a SN la caught ~ 1 week after 
maximum light, or a luminous SN Ic discovered at maximum. There is a weak indication of S II absorption 
at ~ 5375 A which favors classification as a Type la (see Figures 1 and 12), but this alone does not provide 
a secure classification. Note that the if-corrections for a SN la or SN Ic at this redshift (z = 0.43) would 
be nearly identical due to the excellent match of the observed filters (£>45 and to the rest-frame (B 

and V) filters. 

We have reanalyzed the cosmological parameters discarding SN 1996E as a safeguard against the 
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possible contamination of our high-rcdshift sample. We also excluded SN 1997ck which, for lack of a 
definitive spectral classification, is an additional threat to contamination of our sample. With the remaining 
"high-confidence" sample of 14 SNe la we find the statistical likelihood of a positive cosmological constant 
to be 99.8% (3.1 a) from the MLCS method, a modest increase from 99.7% (3.0 a) confidence when SN 
1996E is included. For the template fitting approach, the statistical confidence in a positive cosmological 
constant remains high at >99.9% (4.0 a), the same result as with SN 1996E. We conclude that for this 
sample our results are robust against sample contamination, but the possible contamination of future 
samples remains a concern. Even given existing detector technology, more secure supernova classifications 
can be achieved with greater signal-to-noise ratios for observed spectra, with optimally timed search epochs 
which increase the likelihood of pre-maximum discovery, and with an improved empirical understanding of 
the differences among the spectra of supernova types. 



5.8. Comparisons 

The results reported here are consistent with other reported observations of high-redshift SNe la 
from the High-z Supernova Search Team (Garnavich et al. 1998; Schmidt et al. 1998), and the improved 
statistics of this larger sample reveal the potential influence of a positive cosmological constant. 

These results are inconsistent at the ~ 2a confidence level with those of Perlmutter et al. (1997), 
who found M = 0.94 ± 0.3 (fi A = 0.06) for a flat Universe and M = 0.88 ± 0.64 for fi A = 0. They are 
marginally consistent with those of Perlmutter et al. (1998) who, with the addition of one very high redshift 
SN la (z = 0.83), found fl M = 0.6 ± 0.2 (fi A = 0.4) for a flat Universe and M = 0.2 ± 0.4 for il A = 0. 

Although the experiment reported here is very similar to that performed by Perlmutter et al. (1997, 
1998), there are some differences worth noting. Schmidt et al. (1998), Garnavich et al. (1998), and this 
paper explicitly correct for the effects of extinction evidenced by reddening of the SNe la colors. Not 
correcting for extinction in the nearby and distant sample could affect the cosmological results in either 
direction since we do not know the sign of the difference of the mean extinction. In practice we have found 
few of the high-redshift SNe la to suffer measurable reddening. A number of objects in the nearby sample 
display moderate extinction for which we make individual corrections. We also include the Hubble constant 
as a free parameter in each of our fits to the other cosmological parameters. Treating the nearby sample in 
the same way as the distant sample is a crucial requirement of this work. Our experience observing the 
nearby sample aids our ability to accomplish this goal. 

The statistics of gravitational lenses provide an alternate method for constraining the cosmological 
constant (Turner 1990; Fukugita, Futamase, & Kasai 1990). Although current gravitational lensing 
limits for the cosmological constant in a flat Universe (Q A < 0.66 at 95% confidence; Kochanek 1996) 
are not inconsistent with these results, they are uncomfortably close. Future analysis which seeks to 
limit systematic uncertainties affecting both experiments should yield meaningful comparisons. The 
most incisive independent test may come from measurements of the fluctuation spectrum of the cosmic 
microwave background. While the supernova measurements provide a good constraint on Qm — f2 A , the 
CMB measurements of the angular scale for the first Doppler peak, referring to much earlier epochs, are 
good measures of Om + A (White & Scott 1996). Since these constraints are nearly orthogonal in the 
coordinates of Figure 6 and 7, the region of intersection could be well defined. Ongoing experiments from 
balloons and the South Pole may provide the first clues to the location of where that intersection. 

Our detection of a cosmological constant is not limited by statistical errors but by systematic ones. 
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Further intensive study of SNe la at low (z < 0.1), intermediate (0.1 < z < 0.3), and high (z > 0.3) rcdshifts 
is needed to uncover and quantify lingering systematic uncertainties in this striking result. 

6. Conclusions 

1. We find the luminosity distances to well-observed SNe with 0.16 < z < 0.97 measured by two 
methods to be in excess of the prediction of a low mass-density (£1m « 0.2) Universe by 0.25 to 0.28 mag. 
A cosmological explanation is provided by a positive cosmological constant with 99.7% (3. Oct) to >99.9% 
(4. Oct) confidence using the complete spectroscopic SN la sample and the prior belief that £Im > 0. 

2. The distances to the spectroscopic sample of SNe la measured by two methods are consistent with 
a currently accelerating expansion (q < 0) at the 99.5% (2.8ct) to >99.9% (3.9ct) level for q = - £!a 
using the prior that Qm > 0. 

3. The data favor eternal expansion as the fate of the Universe at the 99.7% (3. Oct) to >99.9% (4. Oct) 
confidence level from the spectroscopic SN la sample and the prior that > 0. 

4. We estimate the dynamical age of the Universe to be 14.2 ±1.5 Gyr including systematic 
uncertainties, but subject to the zeropoint of the current Cepheid distance scale used for the host galaxies 
of three nearby SNe la (Saha et al. 1994, 1997). 

5. These conclusions do not depend on inclusion of SN 1997ck (z=0.97), whose spectroscopic 
classification remains uncertain, nor on which of two light-curve fitting methods is used to determine the 
SN la distances. 

6. The systematic uncertainties presented by grey extinction, sample selection bias, evolution, a local 
void, weak gravitational lensing, and sample contamination currently do not provide a convincing substitute 
for a positive cosmological constant. Further studies are needed to determine the possible influence of any 
remaining systematic uncertainties. 
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A. Appendix: MLCS 

Following the success of Phillips (1993), Riess, Press, & Kirshner (1995) employed a linear estimation 
algorithm (Rybicki & Press 1992) to determine the relationship between the shape of a SN la light curve 
and its peak luminosity. This method was extended (Riess, Press, & Kirshner 1996a) to utilize the SN la 
color curves to quantify the amount of reddening by interstellar extinction. In this Appendix we describe 
further refinements and optimization of the MLCS method for the application to high-redshift SNe la. 

Previously, the MLCS relations were derived from a set similar to the nearby (cz < 2000 km s _1 ) 
sample of Phillips (1993) (Riess, Press, & Kirshner 1995, 1996a). The relative luminosities of this "training 
set" of SNe la were calibrated with independent distance indicators (Tonry 1991; Pierce 1994). The absolute 
SN la luminosities were measured from Cepheid variables populating the host galaxies (Saha et al. 1994, 
1997). Yet at moderate distances, the most reliable distance indicator available in nature is the redshift. 
The recent harvest of SN la samples (Hamuy et al. 1996b; Riess et al. 1998b) with cz > 2500 km s _1 
provides a homogeneous training set of objects for MLCS with well understood relative luminosities. Here 
we employ a set (see Table 10) of B and V light curves with cz > 2500 km s _1 to determine the MLCS 
relations. 

The significant increase in the size of the available training set of SNe la since Riess, Press, & Kirshner 
(1996a) supports an expansion of our description of the MLCS relations. Riess, Press, & Kirshner (1996a) 
described SNe la light curves as a linear family of the peak luminosity: 

m v = M v + RvA + /iy (Al) 

m B -v = Mb-v + Rb-vA + Bb-v (A2) 

where mv,niB-v are the observed light and color curves, A = M v — M v (standard) is the difference in 
maximum luminosity between the fiducial template SN la and any other SN la, Ry and Rb-v are vectors 
of correlation coefficients between A and the light curve shape, fj, v is the apparent distance modulus, and 
Eb-v is the color excess. All symbols in bold denote vectors which are functions of SN la age, with t = 
taken by convention as the epoch of B maximum. 

By adding a second-order term in the expansion, our empirical model becomes 



m v = M v + RvA + Q V A 2 + /xy (A3) 

m B -v = M B -v + Rb-v A + Qb-vA 2 + E B -v (A4) 

where Qv,Qb-v are the correlation coefficients of the quadratic relationship between A 2 and the 
light curve shape. The vectors of coefficients (Ry, Rb-v, Qv, Qb-v) as well as the fiducial templates 
(Mv,Mb-v) are determined from the training set of SNe la listed in Table A. (They can be found at 



http://oir-www.harvard.edu/cfa/oirRcscarch/supcrnova .) The empirical light and color curve families are 



shown in Figure 13. As before, these MLCS relations show that the more luminous SNe la have broader 
light curves and are bluer until day ~ 35, by which time all SNe la have the same color. The primary 
difference from the previous MLCS relations is that near maximum, the color range spanned by the same 
range of SN la luminosities is much reduced. Further, the quadratic MLCS relations reveal that SNe la 
which are brighter or dimmer (than the fiducial value) by equal amounts do not show equal changes in their 
colors. Faint SNe la are far redder than the amount by which luminous SNe la are blue. 

Fitting of this quadratic model (equations A3-A4) to a SN la still requires the determination of 4 
"free" parameters: A, /iy, Eb-v-, and t max . The parameters arc determined by minimizing the expected 



-29 - 



deviations between data and model: 

X 2 = rxC- 1 ^, (A5) 

where 

r x = m x M x - R X A - Q X A 2 - fi x (A6) 

for any band x. Here C is the correlation matrix of the model and the measurements. Correlations of the 
data from the model were determined from the SNe la of Table 10. These correlations result from our still 
imperfect (but improving) description of the light-curve shape behavior. Future expansion of the model will 
reduce these correlations further until they become constraints on the unpredictable, turbulent behavior of 
the SN la atmosphere. Riess, Press, & Kirshner (1996a) quantified the autocorrelation (diagonal matrix 
elements) of the linear model. Here we have determined, in addition, the covariance (off-diagonal matrix 
elements) between two measurements of differing SN la age, passband, or both. These can be found at 



http://oir-www.harvard.edu/cfa/oir/Research/supcrnova. The correlation matrix of the measurements, 



commonly called the "noise," is, as always, provided by the conscientious observer. 

The a priori values for A used to determine the vectors Rv, Rb v, Qv, Qb v, My, and Mb-v 
are the differences between the measured peak magnitudes and those predicted by the SN la host galaxy 
redshift. These values for A must be corrected for the extinction, Ay. Because the values of Ay are not 
known a priori, we use an initial guess derived from the color excess measured from the uniform color range 
of SNe la after day 35 (Riess, Press, & Kirshner 1996a; Lira 1995). 

Initial guesses for A, /Lty, Eb-v, an d t max yield estimates for Rv, Rb v, Qv, Qb v, My, and Mb-v 
by minimizing equation (A5) with respect to the latter. These estimates for Rv, Rb-v, Qv, Qb v, My, 
and Mb-v yield improved estimates of A, [iy , Eb-v, and t max also determined by minimizing equation 
(A5) with respect to the latter. This iterative determination of these vectors and parameters is repeated 
until convergence is reached. Subsequent determination of the parameters A, /j,y, Eb-v, and t max for SNe 
la not listed in Table 10 (such as those reported here) is done using the fixed vectors derived from this 
training process. 

We also employ a refined estimate of the selective absorption to color excess ratio, Rv = Ay / Eb-v , 
which has been calculated explicitly as a function of SN la age from accurate spectrophotometry of SNe la 
(Nugent, Kim, & Perlmutter 1998). This work shows that although Ry is the canonical value of ~ 3.1 for 
SNe la at maximum light or before, over the first 10 days after maximum Ry slowly rises to about 3.4. For 
highly reddened SNe la, this change in Ry over time can appreciably affect the shape of the SN la light 
curve (Leibundgut 1989). 

Lastly wc have refined our a priori understanding of the likelihood for SN la interstellar extinction from 
host galaxies. The previous incarnation of MLCS (Riess, Press, & Kirshner 1996a) employed a "Bayesian 
filter" to combine our measurement of extinction with our prior knowledge of its one-directional effect. In 
addition, it is less probable to observe a very large amount of extinction due to the finite column density of 
a spiral disk as well as a reduced likelihood for detection of SNe with large extinctions. To quantify this a 
priori likelihood for extinction we have adopted the calculations of Hatano, Branch, & Deaton (1997), who 
determined the extinction distribution for SNe la in the bulge and disk of late-type galaxies. The primary 
difference between our previous a priori distribution and the results of Hatano, Branch, & Deaton (1997) 
are that non-trivial quantities of extinction are even less probable than assumed. In particular, Hatano, 
Branch, & Deaton (1998) show that two-thirds of SNe la suffer less than 0.3 to 0.5 mag of extinction, 
which is approximately half the amount of extinction previously assumed. Despite our use of an externally 
derived Bayesian prior for probable SN la extinction, it is important to continue testing that the a posterior 
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cxtinction distribution matches the expected one. A statistically significant departure could imply an 
important deficiency in the SN la luminosity, light-curve shape, and color relations. Specifically, excessively 
blue SNe la such as SN 1994D (E B _ V w —0.10 ± 0.04), if common, would reveal a shortcoming of these 
MLCS relations. However, using the current MLCS relations, the best estimate we can make for such blue 
SNe la is that their extinctions are negligible. If the a priori distributions of Hatano, Branch, & Deaton 
(1997) are not significantly in error, this practice is statistically sensible and does not introduce a distance 
bias. 
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Figure 1: Identification spectra (in i\) of high-redshift SNe la. The spectra obtained for the 10 new 
SNe of the high-redshift sample are shown in the restframe. The data are compared to nearby SN la 
spectra of the same age as determined by the light curves (see Table 1). The spectra the three objects from 
Garnavich et al. (1998) are also displayed. 

Figure 2: Local standard stars in the fields of SNe la. The stars are listed in Table 2 and the locations 
of the stars and SNe are indicated in the figure. The orientation of each field is East to the right and North 
at the top. The width and length of each field is: 96E=4.9', 96H=4.9', 961=4.9', 96J=4.9', 96K=4.9', 
96R=5.0', 96T=4.9', 96U=4.9', 95ao=4.8', 95ap=4.8'. 

Figure 3: Light curves of high-redshift SNe la. B (filled symbols) and V (open symbols) photometry 
in the rest-frame of 10 well-observed SNe la is shown with B increased by 1 mag for ease of view. The lines 
are the empirical MLCS model fits to the data. Supernova age is shown relative to B maximum. 

Figure 4: MLCS SNe la Hubble diagram. The upper panel shows the Hubble diagram for the 
low-redshift and high-redshift SNe la samples with distances measured from the MLCS method (Ricss, 
Press, & Kirshner 1995, 1996a; Appendix of this paper). Overplotted are three cosmologies: "low" and 
"high" Qm with SI a = and the best fit for a flat cosmology, Am = 0.24, Q\ = 0.76. The bottom panel 
shows the difference between data and models with VLm = 0.20, f^A = 0. The open symbol is SN 1997ck 
(z = 0.97) which lacks spectroscopic classification and a color measurement. The average difference between 
the data and the Am = 0.20, Qa = prediction is 0.25 mag. 

Figure 5: Anii5(B) SN la Hubble diagram. The upper panel shows the Hubble diagram for the 
low-redshift and high-redshift SNe la samples with distances measured from the template fitting method 
parameterized by Atoi 5 (£>) (Hamuy et al. 1995, 1996d). Overplotted are three cosmologies: "low" and 
"high" Qm with Q\ = and the best fit for a flat cosmology, CIm = 0.20, Oa = 0.80. The bottom panel 
shows the difference between data and models from the VIm = 0.20, fl\ = prediction. The open symbol 
is SN 1997ck (z = 0.97) which lacks spectroscopic classification and a color measurement. The average 
difference between the data and the fijj = 0.20, Oa = prediction is 0.28 mag. 

Figure 6: Joint confidence intervals for (Dm.^a) from SNe la. The solid contours are results from the 
MLCS method applied to well-observed SNe la light curves together with the snapshot method (Riess et al. 
1998a) applied to incomplete SNe la light curves. The dotted contours are for the same objects excluding 
the unclassified SN 1997ck (z — 0.97). Regions representing specific cosmological scenarios are illustrated. 
Contours are closed by their intersection with the line Q.m = 0. 

Figure 7: Joint confidence intervals for (Om^a) from SNe la. The solid contours are results from the 
template fitting method applied to well-observed SNe la light curves together with the snapshot method 
(Ricss et al. 1998a) applied to incomplete SNe la light curves. The dotted contours are for the same objects 
excluding the unclassified SN 1997ck (z = 0.97). Regions representing specific cosmological scenarios are 
illustrated. Contours are closed by their intersection with the line CIm = 0. 

Figure 8: PDF for the dynamical age of the Universe from SNe la (equation 19). The PDF for the 
dynamical age derived from the PDFs for H , QmMa is shown for the two different distance methods 
without the unclassified SN 1997ck. A naive average (see §4.2) yields an estimate of 14.2+^g Gyr, not 
including the systematic uncertainties in the Ccphcid distance scale. 

Figure 9: Lines of constant dynamical age in Gyr in the (Um^a) plane. Comparing these lines with 
the error ellipses in Figures 5 and 7 reveals the leverage this experiment has on measuring the dynamical 
age. This plot assumes H = 65 km s _1 Mpc -1 as determined from nearby SNe la and is subject to the 
zeropoint of the Cepheid distance scale. 

Figure 10: Distributions of MLCS light curve shape parameters, A, and template fitting parameters, 
Ami 5 (B), for the high and low-redshift samples of SNe la. Positive values for A and Amis(B) > 1.1 
correspond to intrinsically dim SNe la, negative values for A and Atoi 5 (£>) < 1.1 correspond to luminous 
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SNe la. Histograms of the low-redshift (solid line) and high-redshift (dotted line) light curve shape 
parameters are mutually consistent with no indication that these samples arc drawn from different 
populations of SNe la. Filled and open circles show the distribution of log(cz) for the low and high-redshift 
samples, respectively. 

Figure 11: Spectral comparison (in f\) of SN 1998ai (z — 0.49) with low-redshift (z < 0.1) SNe la 
at a similar age. Within the narrow range of SN la spectral features, SN 1998ai is indistinguishable from 
the low-redshift SNe la. The spectra from top to bottom are SN 1992A, SN 1994B, SN 1995E, SN 1998ai, 
and SN 1989B ~ 5 days before maximum light. The spectra of the low-redshift SNe la were resampled and 
convolved with Gaussian noise to match the quality of the spectrum of SN 1998ai. 

Figure 12: Comparison of the spectral and photometric observations of SN 1996E to those of Type 
la and Type Ic supernovae. The low signal-to-noise ratio of the spectrum of SN 1996E and the absence of 
data blueward of 4500 A makes it difficult to distinguish between a Type la and Ic classification. The light 
and color curves of SN 1996E are also consistent with either supernova type. The spectrum was taken six 
days (rest-frame) after the first photometric observation. 

Figure 13: MLCS empirical SN la light curve families in M B , M v , and (B - V)q. The derived light 
curves are given as a function of the luminosity difference, A, between the peak visual luminosity of a SN la 
and a fiducial (A = 0) SN la. Properties of the SN la families are indicated in the figure and the Appendix. 
The light and color curves of SN 1995ac (open symbols) and SN 1996X (filled symbols) are overplotted as 
examples of luminous and dim SN la, respectively. 



This figure "ariess.figl.gif" is available in "gif" format from: 
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Table 8: Cosmological Results 













no constraint* 






Qtot = 1 


Q A = 


fi M = 0.2 


Method (high-z SNe) 


Ho 








to 


p(Sl A > 0) 


p(qo < 0) 


go 








T MLCS±Snap.(15) 








1.19 




99.7%(3.0cr) 


99.5%(2.8<r) 


-0.98±0.40 


0.28±0.10 


-0.34±0.21 


0.65 ±0.22 


T AMi 5 ±Snap.(15) 








1.03 




>99.9%(4.0cr) 


>99.9%(3.9<r) 


-1.34±0.40 


0.17±0.09 


-0.48±0.19 


0.84 ±0.18 


MLCS+Snap.+97ck(16) 




0.24l° f 4 


72+° 72 

"• 1 ^-0.48 


1.17 




99.5%(2.8<r) 


99.3%(2.7<r) 


-0.75±0.32 


0.24±0.10 


-0.35±0.18 


0.66 ±0.21 


AMi5+Snap.+97ck(16) 




0.80t° ;t° 8 


1 -° {i —0.70 


1.04 




>99.9%(3.9<r) 


>99.9%(3.8<r) 


-1.14±0.30 


0.21±0.09 


-0.41±0.17 


0.80 ±0.19 


MLCS(9) 


65.2±1.3 { 






1.19 


13 fi +1 


99.6%(2.9<r) 


99.4%(2.4<r) 


-0.92±0.42 


0.28±0.10 


-0.38±0.22 


0.68 ±0.24 


AAfi 5 (9) 


63.8±1.3 { 






1.05 


14 8 +1 

1 ^-°-0.8 


>99.9%(3.9<r) 


>99.9%(3.8<r) 


-1.38±0.46 


0.16±0.09 


-0.52±0.20 


0.88 ±0.19 


MLCS+97ck(10) 


65.2±1.3 { 


00+° 60 

"•""-0.00 


48+° 72 

"•^°-0.24 


1.17 


14.2li 3 


99.5%(2.8<r) 


99.3%(2.7<r) 


-0.74±0.32 


0.24±0.10 


-0.38±0.19 


0.68 ±0.22 


AMi5+97ck(10) 


63.7±1.3 { 


72+° 44 

"• 1 ^-0.56 


1 48+° 56 


1.04 


J-O- J--0.9 


>99.9%(3.8<r) 


> 99.9%(3.7<r) 


-1.11±0.32 


0.20±0.09 


-0.44±0.18 


0.84 ±0.20 


Snap. (6) 


63.4±2.7 { 






1.30 




89.1%(1.6<r) 


78.9%(1.3cr) 


-0.70±0.80 


0.40±0.50 


0.06 ±0.70 


0.44 ±0.60 



*fi M > 

T Complete set of spectroscopic SNe la. 

*This uncertainty reflects only the statistical error from the variance of SNe la in the Hubble flow. 
It does not include any contribution from the (much larger) SN la absolute magnitude error. 
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Table 1. High-z Supernova Spectroscopy 



SN 


UT Date 


Telescope 


Spectral Range (nm) 


Redshift 


Comparison 


1995ao 


95 Nov 23 


Keck-I 


510-1000 


0.24 a 


1996X(-4) 


1995ap 


95 Nov 23 


Keck-I 


510-1000 


0.30 b 


1996X(-4) 


1996E 


96 Feb 23 


ES03.6m 


600-990 


0.43 a 


1989B(+9) 


1996H 


96 Feb 23 


ES03.6m 


600-990 


0.62 a 


1996X(+5) 


19961 


96 Feb 23 


ES03.6m 


600-990 


0.57 b 


1996X(+5) 


1996J 


96 Feb 23 


ES03.6m 


600-990 


0.30 a 


1995D(+0) 


1996K 


96 Feb 23 


ES03.6m 


600-990 


0.38 b 


1995D(+0) 


1996R 


96 Mar 18 


MMT 


400-900 


0.16 a 


1989B(+12) 


1996T 


96 Mar 18 


MMT 


400-900 


0.24 a 


1996X(-4) 


1996U 


96 Mar 18 


MMT 


400-900 


0.43 a 


1995D(+0) 


1997ce 


97 May 04 


Keck-II 


570-940 


0.44 b 


1995D(+0) 


1997cj 


97 May 02 


MMT 


400-900 


0.50 a 




1997cj 


97 May 04 


Keck-II 


570-940 


0.50 b 


1995D(+0) 


1997ck 


97 May 04 


Keck-II 


570-940 


0.97 a 





a Derived from emission lines in host galaxy. 
b Derived from broad features in SN spectrum. 

c Supernova and its age (relative to B maximum) used for comparison spectrum in Figure 1. 
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Table 2. SN la Imaging 



JD a 


UT Date 


B45 


V45 


B35 


V35 


Telescope 








SN 1996E 








127.6 


1996 Feb 14 


22.30(0.09) 








CTIO 4m 


128.6 


1996 Feb 15 


22.27(0.04) 


21.86(0.08) 






CTIO 4m 


132.1 


1996 Feb 19 


22.46^(0.11) 








ESO NTT 


134.6 


1996 Feb 21 


22.66(0.10) 


21.99(0.26) 






CTIO 4m 


135.5 


1996 Feb 22 


22.68(0.13) 


22.09(0.06) 






CTIO 4m 


138.7 


1996 Feb 25 


23.04(0.12) 


22.29(0.15) 






ESO 1.5m 


139.6 


1996 Feb 26 


22.89(0.15) 


22.72(0.33) 






ESO 1.5m 


157.6 


1996 Mar 15 


24.32(0.18) 


23.51(0.77) 






CTIO 4m 


163.7 


1996 Mar 21 




22.87(0.50) 
SN 1996H 






WIYN 


127.6 


1996 Feb 14 


22.78(0.13) 








CTIO 4m 


128.6 


1996 Feb 15 


22.81(0.06) 


22.25(0.14) 






CTIO 4m 


132.1 


1996 Feb 19 


22.71^(0.29) 


22.40 7 (0.37) 






ESO NTT 


134.6 


1996 Feb 21 


22.85(0.08) 


22.48(0.19) 






CTIO 4m 


135.5 


1996 Feb 22 


22.83(0.18) 


22.28(0.10) 






ESO 3.6m 


136.6 


1996 Feb 23 


22.84(0.13) 








ESO 3.6m 


138.7 


1996 Feb 25 


22.85(0.09) 


22.58(0.15) 






ESO 1.5m 


139.6 


1996 Feb 26 


22.88(0.15) 


22.52(0.25) 






ESO 1.5m 


140.6 


1996 Feb 27 


22.96(0.16) 


23.10(0.10) 






ESO 1.5m 


141.6 


1996 Feb 28 


23.05(0.08) 








WIYN 


142.6 


1996 Feb 29 


23.21(0.20) 


22.69(0.16) 






WIYN 


157.6 


1996 Mar 15 


23.98(0.22) 


23.18(0.28) 






CTIO 4m 


161.6 


1996 Mar 19 


24.16(0.22) 








CTIO 4m 


164.6 


1996 Mar 22 




24.01(0.30) 






WIYN 
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Table 2 — Continued 



JD a 


UT Date 


B45 


1/45 


B35 


V35 


Telescope 








SN 19961 








128.6 


1996 Feb 15 


22.77(0.05) 








CTIO 4m 


132.1 


1996 Feb 19 


22.95(0.22) 


22.30(0.22) 






ESO NTT 


134.6 


1996 Feb 21 


22.95(0.05) 


22.65(0.15) 






CTIO 4m 


135.5 


1996 Feb 22 


22.92(0.05) 


22.64(0.20) 






ESO 3.6m 


136.6 


1996 Feb 23 


22.88(0.09) 


22.74(0.28) 






ESO 3.6m 


138.7 


1996 Feb 25 


23.12(0.13) 


22.86(0.17) 






ESO 1.5m 


140.6 


1996 Feb 27 


23.64(0.36) 


22.67(0.36) 






ESO 1.5m 


142.6 


1996 Feb 29 


23.48(0.10) 


23.06(0.22) 






WIYN 


157.6 


1996 Mar 15 


24.83(0.17) 


23.66(0.30) 






CTIO 4m 


161.6 


1996 Mar 19 


24.70(0.31) 


SN 1996J 






CTIO 4m 


127.6 


1996 Feb 14 


22.01(0.02) 








CTIO 4m 


128.6 


1996 Feb 15 


21.95(0.03) 


21.95(0.07) 






CTIO 4m 


134.6 


1996 Feb 21 


21.57(0.03) 


21.59(0.05) 






CTIO 4m 


135.6 


1996 Feb 22 


21.62(0.04) 


21.61(0.04) 


21.84(0.03) 


21.46(0.06) 


ESO 3.6m 


135.6 


1996 Feb 22 






21.89(0.04) 


21.47(0.02) 


ESO 3.6m 


139.7 


1996 Feb 26 


21.63(0.04) 


21.46(0.07) 






ESO 1.5m 


140.7 


1996 Feb 27 






21.90(0.07) 


21.77(0.05) 


ESO 1.5m 


157.6 


1996 Mar 15 


22.77(0.05) 


22.06(0.12) 


23.69(0.07) 


21.83(0.04) 


CTIO 4m 


161.8 


1996 Marl9 






24.34(0.19) 


22.05(0.05) 


CTIO 4m 


166.6 


1996 Mar 24 








22.76(0.07) 


CTIO 1.5m 
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Table 2 — Continued 



JD a 


UT Date 


B45 


1/45 


B35 


V35 


Telescope 








SN 1996K 








128.5 


1996 Feb 15 


23.74(0.04) 








CTIO 4m 


135.5 


1996 Feb 22 


22.49(0.07) 








ESO 3.6m 


135.5 


1996 Feb 22 


22.52(0.07) 








ESO 3.6m 


135.7 


1996 Feb 22 


22.56(0.03) 


22.48(0.06) 






ESO 3.6m 


136.6 


1996 Feb 23 


22.48(0.05) 


22.26(0.16) 






ESO 1.5m 


138.6 


1996 Feb 25 


22.15(0.10) 


22.47(0.11) 






ESO 1.5m 


138.7 


1996 Feb 25 


22.18(0.07) 








ESO 1.5m 


139.6 


1996 Feb 26 


22.37(0.05) 


22.42(0.13) 






ESO 1.5m 


140.8 


1996 Feb 27 






22.23(0.10) 


22.06(0.11) 


ESO 1.5m 


157.5 


1996 Mar 15 


22.83(0.07) 




22.93(0.12) 


22.61(0.19) 


CTIO 4m 


157.5 


1996 Mar 15 


22.81(0.09) 




22.86(0.10) 


22.45(0.10) 


CTIO 4m 


161.7 


1996 Mar 19 


23.20(0.16) 


22.45(0.13) 


23.17(0.17) 


22.69(0.15) 


CTIO 4m 


162.6 


1996 Mar 20 


23.17(0.06) 


22.79(0.12) 






WIYN 


165.6 


1996 Mar 23 






23.58(0.16) 


23.17(0.14) 


CTIO l.m 


168.5 


1996 Mar 26 








23.20(0.19) 


CTIO l.m 


169.7 


1996 Mar 27 


24.05(0.26) 


SN 1996R 


24.42(0.25) 




MDM 


157.7 


1996 Mar 15 


20.48(0.01) 








CTIO 4m 


158.7 


1996 Mar 16 


20.59(0.03) 


20.70(0.03) 






CTIO 4m 


167.7 


1996 Mar 25 








21.62 v '(0.04) 


CTIO 1.5m 


191.7 


1996 Apr 18 


22.41(0.09) 


SN 1996T 






ESO 1.5m 


161.7 


1996 Mar 19 


20.83^(0.03) 




20.86 v '(0.02) 




CTIO 4m 


167.6 


1996 Mar 25 


20.95^(0.04) 




20.96 v '(0.03) 




CTIO 1.5m 


191.7 


1996 Apr 18 






22.37 F (0.17) 




ESO 1.5m 


212.6 


1996 May 9 


22.52^(0.08) 




22.99 F (0.31) 




WIYN 
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Table 2 — Continued 



JD a 


UT Date 


B45 


1/45 


B35 


V35 


Telescope 








SN 1996U 








158.7 


1996 Mar 16 


22.16(0.04) 








CTIO 4m 


160.7 


1996 Mar 18 


22.00(0.11) 


22.03(0.18) 






MDM 


161.7 


1996 Mar 19 


22.04(0.05) 


22.23(0.26) 






CTIO 4m 


165.7 


1996 Mar 23 




22.35(0.28) 






CTIO 1.5m 


167.7 


1996 Mar 25 


22.19(0.10) 








CTIO 1.5m 


186.7 


1996 Apr 13 


23.33^(0.17) 


22.64 7 (0.28) 






LCO 


188.7 


1996 Apr 15 


23.51(0.17) 


22.96(0.36) 

SN 1995ao 






WIYN 


39.6 


1995 Nov 18 


21.42(0.05) 








CTIO 4m 


46.6 


1995 Nov 25 


21.30(0.03) 


21.10(0.13) 






WIYN 


51.6 


1995 Nov 30 


21.24(0.05) 


SN 1995ap 


21.52(0.05) 


21.12(0.03) 


CTIO 4m 


39.6 


1995 Nov 18 


22.41(0.14) 








CTIO 4m 


46.6 


1995 Nov 25 


21.13(0.08) 


21.40(0.10) 






WIYN 


48.6 


1995 Nov 27 


21.04(0.11) 








WIYN 


51.6 


1995 Nov 30 


21.04(0.11) 




21.65(0.09) 


20.92(0.07) 


CTIO 4m 



Note. — a actually JD-2450000. 
Uncertainties in magnitudes are listed in parentheses. 
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Table 3. SN la Field Local Standard Stars 



Star B45 V45 B35 V35 



SN 1996E 

1 20.84(0.02) 20.71(0.02) 

2 20.07(0.03) 18.69(0.03) 

3 19.60(0.03) 19.22(0.03) 

4 19.76(0.03) 18.35(0.03) 

5 19.16(0.03) 18.29(0.03) 

6 20.85(0.02) 20.52(0.02) 

SN 1996H 

1 18.16(0.02) 17.84(0.02) 

2 19.96(0.02) 18.50(0.02) 

3 21.13(0.02) 19.41(0.02) 

4 20.76(0.02) 19.21(0.02) 

5 19.62(0.02) 19.23(0.02) 

6 20.02(0.02) 19.69(0.02) 

SN 19961 

1 19.59(0.02) 18.67(0.02) 

2 22.35(0.02) 20.72(0.02) 

3 20.62(0.02) 18.93(0.02) 

4 20.22(0.02) 18.97(0.02) 

5 17.46(0.02) 17.18(0.02) 

6 18.02(0.02) 17.55(0.02) 
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Table 3 — Continued 



Star 


B45 


1/45 


B35 


V35 






SN 1996J 






1 


18.59(0.02) 


17.38(0.02) 


19.09(0.02) 


17.85(0.02) 


2 


20.27(0.02) 


19.49(0.02) 


20.74(0.02) 


19.78(0.02) 


3 


20.20(0.02) 


19.45(0.02) 


20.70(0.02) 


19.79(0.02) 


4 


19.63(0.02) 


18.67(0.02) 


20.09(0.02) 


19.06(0.02) 


5 


21.12(0.02) 


19.63(0.02) 


21.69(0.02) 


20.20(0.02) 


6 


20.27(0.02) 


20.00(0.02) 


20.57(0.02) 


20.06(0.02) 






SN 1996K 






1 


19.06(0.02) 


18.81(0.02) 


19.22(0.02) 


18.88(0.02) 


2 


19.76(0.03) 


19.43(0.03) 


19.94(0.03) 


19.53(0.03) 


3 


19.41(0.03) 


18.17(0.02) 


19.90(0.03) 


18.62(0.02) 


4 


19.84(0.03) 


18.64(0.02) 


20.28(0.03) 


19.06(0.02) 


5 


19.30(0.02) 


17.70(0.02) 


19.84(0.02) 


18.25(0.02) 


6 


19.04(0.02) 


18.06(0.02) 


19.45(0.02) 


18.40(0.02) 


7 


18.05(0.02) 


17.17(0.02) 


18.47(0.02) 


17.49(0.02) 






SN 1996R 






1 


17.29(0.03) 


16.61(0.02) 


18.01 ,/ (0.03) 




2 


18.15(0.03) 


17.78(0.03) 


is^s^o.os) 




3 


19.05(0.03) 


18.67(0.03) 


lg.si^o.os) 




4 


19.20(0.03) 


18.22(0.03) 


20.02 v '(0.03) 




5 


18.06(0.03) 


17.64(0.03) 


18.54 ,/ (0.03) 
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Table 3 — Continued 



Star 


B45 


V45 


B35 


V35 






SN 1996T 






1 


18.29(0.02) v 


18.01(0.02)^ 






2 


19.77(0.02) v 


18.57(0.02)^ 






3 


20.43(0.02) ,/ 


19.31(0.02) H 






4 


21.28(0.02) ,/ 


20.57(0.02)^ 






5 


21.28(0.02) ,/ 


20.27(0.02)^ 






6 


21.34(0.02) ,/ 


20.37(0.02)^ 

SN 1995ao 






1 


20.36(0.03) 


20.15(0.03) 


20.59(0.03) 


20.19(0.03) 


2 






17.89(0.03) 


17.50(0.03) 


3 


20.09(0.03) 


19.50(0.03) 


20.50(0.03) 


19.79(0.03) 


4 


20.10(0.03) 


19.75(0.03) 


20.39(0.03) 


19.86(0.03) 


5 


16.37(0.03) 


15.47(0.03) 


16.62(0.03) 


15.73(0.03) 


6 




SN 1995ap 


17.32(0.03) 


16.81(0.03) 


1 


19.49(0.03) 


18.21(0.03) 


20.28(0.03) 


18.69(0.02) 


2 


19.19(0.03) 


18.76(0.03) 


19.54(0.03) 


18.88(0.02) 


3 


18.97(0.03) 


18.24(0.03) 


19.43(0.03) 


18.47(0.02) 


4 


19.67(0.03) 


18.61(0.03) 


20.31(0.02) 


18.98(0.02) 


5 


20.51(0.03) 


19.44(0.03) 


21.16(0.02) 


19.81(0.02) 


6 


20.90(0.03) 


20.31(0.03) 


21.53(0.02) 


20.50(0.02) 



Note. — Uncertainties in magnitudes are listed in parentheses. 



Table 4. SNe la Light Curves Corrected to the Rest-Frame 



JD a 



SN 1996E 



127.6 


23.04(0.09) 




-0.74 




128.6 


23.02(0.05) 


22.72(0.08) 


-0.74 


-0.86 


132.1 


23.23(0.11) 




-0.77 




134.6 


23.39(0.10) 


22.84(0.26) 


-0.73 


-0.85 


135.5 


23.41(0.13) 


22.95(0.06) 


-0.74 


-0.85 


138.7 


23.80(0.12) 


23.14(0.15) 


-0.76 


-0.84 


139.6 


23.65(0.15) 


23.57(0.33) 


-0.76 


-0.84 


157.6 


25.14(0.18) 


24.42(0.77) 


-0.82 


-0.91 


163.7 




23.78(0.50) 
olN 199611 




-0.91 


127.6 


23.32(0.13) 




-0.54 




128.6 


o o on/ n nn\ 

23.39(0.09) 


23.42(0.14) 


-0.58 


-1.17 


132.1 


23.27(0.30) 


23.56(0.37) 


-0.56 


-1.15 


134.6 


23.48(0.11) 


23.58(0.19) 


-0.63 


-1.10 


135.5 


23.29(0.18) 


23.40(0.10) 


-0.47 


-1.12 


136.6 


23.29(0.14 ) 




-0.45 




138.7 


23.47(0.10) 


23.64(0.15) 


-0.62 


-1.06 


139.6 


23.55(0.18) 


23.56(0.25) 


-0.67 


-1.04 


140.6 


23.58(0.18) 


24.13(0.11) 


-0.62 


-1.03 


141.6 


23.62(0.12) 




-0.56 




142.6 


23.74(0.21) 


23.72(0.16) 


-0.53 


-1.04 


157.6 


O A A A i C\ O O \ 

24.44(0.22) 


23.86(0.28) 


-0.47 


-0.69 


161.6 


24.60(0.22) 




-0.44 




164.6 




24.62(0.30) 
oiN 19961 




-0.61 


128.6 


23.45(0.08) 




-0.68 




132.1 


23.62(0.22) 


23.25(0.22) 


-0.67 


-0.96 


134.6 


23.57(0.06) 


23.66(0.15) 


-0.62 


-1.02 


135.5 


23.54(0.06 ) 


23.65(0.20 ) 


-0.62 


-1.01 


136.6 


23.58(0.10) 


23.74(0.28) 


-0.70 


-1.00 


138.7 


23.81(0.14) 


23.82(0.17) 


-0.69 


-0.97 


140.6 


24.28(0.36 ) 


23.66(0.36 ) 


-0.64 


-1 .00 


142.6 


24.13(0.10 ) 


24.02(0.22 ) 


-0.65 


-0.97 


157.6 


25.38(0.17) 


24.39(0.30) 


-0.55 


-0.73 


161.6 


25.25(0.31) 


SN 1996J 


-0.55 




127.6 


22.58(0.03) 




-0.57 




128.6 


22.52(0.04) 


22.72(0.07) 


-0.57 


-0.78 


134.6 


22.22(0.03) 


22.27(0.06) 


-0.64 


-0.68 


135.6 


22.34(0.03) 


22.08(0.06) 


-0.50 


-0.62 


135.6 


22.39(0.04) 


22.08(0.06) 


-0.50 


-0.62 


135.6 


22.27(0.05) 


22.08(0.06) 


-0.65 


-0.62 


139.7 


22.31(0.06) 


22.15(0.09) 


-0.68 


-0.69 


140.7 


22.41(0.07) 


22.39(0.05) 


-0.51 


-0.62 


157.6 


24.17(0.15) 


22.52(0.04) 


-0.62 


-0.69 


157.6 


23.86(0.06) 


22.52(0.04) 


-1.09 


-0.69 


161.8 




22.76(0.05) 




-0.71 


168.7 




23.48(0.07) 

SN 1996K 




-0.73 


128.5 


24.24(0.05) 




-0.50 




135.5 


23.02(0.07) 




-0.53 
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Table 4 — Continued 



JD a 


B 


V 


K B 


K v 


135.5 


23.04(0.07) 




-0.53 




135.7 


23.09(0.03) 


23.19(0.06) 


-0.53 


-0.71 


136.6 


23.01(0.05) 


22.98(0.16) 


-0.53 


-0.72 


138.6 


22.68(0.10) 


23.20(0.11) 


-0.53 


-0.72 


138.6 


22.71(0.07) 


23.20(0.11) 


-0.53 


-0.72 


139.6 


22.90(0.05) 


23.15(0.13) 


-0.54 


-0.73 


140.8 


22.76(0.10) 


22.77(0.11) 


-0.53 


-0.71 


157.5 


23.43(0.12) 


23.29(0.19) 


-0.50 


-0.67 


157.5 


23.35(0.10) 


23.29(0.19) 


-0.50 


-0.67 


157.5 


23.61(0.07) 


23.29(0.19) 


-0.79 


-0.67 


157.5 


23.59(0.09) 


23.29(0.19) 


-0.79 


-0.67 


161.7 


23.96(0.16) 


23.23(0.13) 


-0.77 


-0.78 


161.7 


23.68(0.17) 


23.23(0.13) 


-0.51 


-0.78 


162.6 


23.95(0.07) 


23.57(0.12) 


-0.79 


-0.78 


165.6 


24.09(0.16) 


23.80(0.14) 


-0.50 


-0.62 


168.5 




23.73(0.19) 




-0.63 


169.6 


24.92(0.25) 




-0.50 




169.7 


24.91(0.26) 


SN 1996R 


-0.86 




157.7 




20.81(0.02) 




-0.33 


158.7 




20.92(0.03) 




-0.33 


167.7 


22.24(0.03) 




-0.63 




191.7 




22.76(0.09) 

SN 1996T 




-0.35 


161.7 


21.24(0.02) 


21.27(0.03) 


-0.38 


-0.44 


167.6 


21.34(0.03) 


21.35(0.04) 


-0.38 


-0.40 


191.7 


22.73(0.20) 




-0.37 




212.6 


23.35(0.35) 


22.81(0.09) 

SN 1996U 


-0.36 


-0.29 


158.7 


22.89(0.05) 




-0.73 




160.7 


22.73(0.11) 


22.88(0.18) 


-0.73 


-0.85 


161.7 


22.78(0.05) 


23.09(0.26) 


-0.74 


-0.85 


165.7 




23.21(0.28) 




-0.86 


167.7 


22.94(0.10) 




-0.75 




186.7 


24.23(0.17) 


23.49(0.28) 


-0.90 


-0.86 


188.7 


24.34(0.17) 


23.85(0.36) 

SN 1995ao 


-0.83 


-0.89 


46.6 


21.85^(0.13) 




-0.75 




51.6 


21.95(0.05) 


21.70(0.03) 

SN 1995ap 


-0.43 


-0.58 


39.6 




22.85(0.14) 




-0.44 


46.6 


21.96^(0.10) 


21.57(0.08) 


-0.56 


-0.45 


48.6 




21.49(0.11) 




-0.45 


51.6 


21.84(0.09) 


21.40(0.08) 


-0.20 


-0.89 



Note. — a actually JD-2450000. 
Uncertainties in magnitudes are listed in parentheses. 
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Table 5. High-z MLCS SN la Light Curve Parameters 



SN 


z 


max 
" l B 


max 

Illy 


A 


A B 


MO (<^Mo) 


1996E 


0.43 


22.81(0.21) 


22.72(0.23) 


-0.08(0.19) 


0.31 


41.74(0.28) 


1996H 


0.62 


23.23(0.19) 


23.56(0.18) 


-0.42(0.16) 


0.00 


42.98(0.17) 


19961 


0.57 


23.35(0.28) 


23.59(0.26) 


-0.06(0.26) 


0.00 


42.76(0.19) 


1996J 


0.30 


22.23(0.12) 


22.21(0.11) 


-0.22(0.10) 


0.24 


41.38(0.24) 


1996K 


0.38 


22.64(0.12) 


22.84(0.14) 


0.29(0.06) 


0.00 


41.63(0.20) 


1996U 


0.43 


22.78(0.22) 


22.98(0.30) 


-0.52(0.29) 


0.00 


42.55(0.25) 


1997ce 


0.44 


22.85(0.09) 


22.95(0.09) 


0.07(0.08) 


0.00 


41.95(0.17) 


1997cj 


0.50 


23.19(0.11) 


23.29(0.12) 


-0.04(0.11) 


0.00 


42.40(0.17) 


1997ck 


0.97 


24.78(0.25) 




-0.19(0.23) 




44.39(0.30) 


1995K 


0.48 


22.91(0.13) 


23.08(0.20) 


-0.33(0.26) 


0.00 


42.45(0.17) 



Note. — Uncertainties in magnitudes are listed in parentheses. 



Table 6. High-z Template Fitting SN la Light Curve Parameters 



SN 


z 


max 
" l B 


max 

Illy 


AMi 5 (B) 


A B 


MO (<^Mo) 


1996E 


0.43 


22.72(0.19) 


22.60(0.12) 


1.18(0.13) 


0.10 


42.03(0.22) 


1996H 


0.62 


23.31(0.06) 


23.57(0.06) 


0.87(0.05) 


0.00 


43.01(0.15) 


19961 


0.57 


23.42(0.08) 


23.61(0.08) 


1.39(0.17) 


0.00 


42.83(0.21) 


1996J 


0.30 


22.28(0.05) 


22.06(0.05) 


1.27(0.27) 


0.64 


40.99(0.25) 


1996K 


0.38 


22.80(0.05) 


22.86(0.08) 


1.31(0.14) 


0.00 


42.21(0.18) 


1996U 


0.43 


22.77(0.05) 


22.96(0.11) 


1.18(0.10) 


0.00 


42.34(0.17) 


1997ce 


0.44 


22.83(0.05) 


22.92(0.05) 


1.30(0.06) 


0.00 


42.26(0.16) 


1997cj 


0.50 


23.29(0.05) 


23.29(0.05) 


1.16(0.03) 


0.09 


42.70(0.16) 


1997ck 


0.97 


24.78(0.16) 




1.00(0.17) 




44.30(0.19) 


1995K 


0.48 


22.92(0.08) 


23.07(0.07) 


1.16(0.18) 


0.00 


42.49(0.17) 



Note. — Uncertainties in magnitudes are listed in parentheses. 
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Table 7. High-z SN la Snapshot Parameters 



SN 


z 


t spec 


A 


A v 


MO (<^Mo) 


1995ao 


0.30 


-2.8 


0.35 


0.00 


40.74(0.60) 


1995ap 


0.23 


-2.9 


0.69 


0.00 


40.33(0.46) 


1996R 


0.16 


8.6 


0.28 


0.10 


39.08(0.40) 


1996T 


0.24 


-4.5 


-0.12 


0.10 


40.68(0.43) 


19971* 


0.17 


0.1 


-0.39 


0.00 


39.95(0.24) 


1997ap* 


0.83 


-2.0 


0.00 


0.00 


43.67(0.35) 



Note. — * as given in Riess et al. 1998a (see also Perlmutter et al. 1998) 



Table 8. Cosmological Results: See Tables8.tex 



Method 


Ho 








to 


p(fi A > 0) 









Table 9. Nearby SN la Snapshot Parameters* 



SN 


log cz 


A 


Ay fit 


M0 (<^Mo) 


1994U 


3.111 


0.03 


0.70 


31.72(0.10) 


1997bp 


3.363 


-0.26 


0.62 


32.81(0.10) 


1996V 


3.870 


0.26 


0.00 


35.35(0.17) 


1994C 


4.189 


0.81 


0.00 


36.72(0.15) 


1995M 


4.202 


-0.15 


0.46 


37.12(0.15) 


1995ae 


4.308 


0.38 


0.00 


37.58(0.21) 


1994B 


4.431 


-0.02 


0.38 


38.51(0.10) 



Note. — 



* as given in Riess et al. 1998a 
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Table 10: Nearby MLCS and Template Fitting SN la Parameters 
MLCS template fitting 



SN 


log cz 


A 


A B 


Mo 


(*) 






Ami 5 (B) 


A B 




(*) 






Host Type* 


1992bo 


3.734 


0.31 


0.00 


34 


72 


;o 


16) 


1.59 


0.01 


34 


88 


;o 


21) 


E 


1992bc 


3.779 


-0.50 


0.00 


34 


87 


;o 


ii) 


0.88 


0.00 


34 


77 


;o 


15) 


L 


1992aq 


4.481 


0.05 


0.00 


38 


41 


;o 


15) 


1.12 


0.18 


38 


33 


;o 


23) 


L 


1992ae 


4.350 


-0.05 


0.00 


37 


80 


;o 


17) 


1.21 


0.28 


37 


77 


;o 


19) 


E 


1992P 


3.896 


-0.19 


0.00 


35 


76 


;o 


13) 


0.94 


0.14 


35 


59 


;o 


16) 


L 


1990af 


4.178 


0.09 


0.18 


36 


53 


;o 


15) 


1.66 


0.16 


36 


67 


;o 


25) 


L 


1994M 


3.859 


0.04 


0.08 


35 


39 


;o 


18) 


1.47 


0.06 


35 


49 


;o 


20) 


E 


1994S 


3.685 


-0.44 


0.00 


34 


27 


;o 


12) 


1.02 


0.03 


34 


34 


;o 


14) 


L 


1994T 


4.030 


0.11 


0.22 


36 


19 


;o 


21) 


1.35 


0.19 


36 


50 


;o 


20) 


L 


1995D 


3.398 


-0.42 


0.00 


33 


01 


;o 


13) 


0.96 


0.23 


32 


79 


;o 


16) 


E 


1995E 


3.547 


-0.61 


2.67 


33 


60 


;o 


17) 


1.03 


2.47 


33 


73 


;o 


17) 


L 


1995ac 


4.166 


-0.47 


0.00 


36 


85 


;o 


13) 


0.99 


0.22 


36 


60 


;o 


16) 


L 


1995ak 


3.820 


0.15 


0.00 


35 


15 


;o 


16) 


1.28 


0.06 


35 


43 


;o 


18) 


L 


1995bd 


3.679 


-0.29 


2.52 


34 


15 


;o 


19) 


0.87 


2.67 


34 


00 


;o 


18) 


L 


1996C 


3.924 


-0.07 


0.24 


35 


98 


;o 


20) 


0.97 


0.42 


35 


82 


;o 


20) 


L 


1996ab 


4.572 


-0.13 


0.00 


39 


01 


;o 


13) 


1.10 


0.09 


39 


10 


;o 


17) 


L 


1992ag 


3.891 


-0.50 


0.77 


35 


37 


;o 


23) 


1.12 


0.09 


35 


53 


;o 


20) 


L 


1992al 


3.625 


-0.35 


0.00 


33 


92 


;o 


11) 


1.13 


0.00 


34 


13 


;o 


14) 


L 


1992bg 


4.024 


-0.06 


0.50 


36 


26 


;o 


21) 


1.14 


0.44 


36 


49 


;o 


21) 


L 


1992bh 


4.130 


-0.16 


0.28 


36 


91 


;o 


17) 


1.05 


0.34 


36 


87 


;o 


17) 


L 


1992bl 


4.111 


-0.06 


0.00 


36 


26 


;o 


15) 


1.50 


0.00 


36 


53 


;o 


20) 


L 


1992bp 


4.379 


-0.26 


0.04 


37 


65 


;o 


13) 


1.27 


0.03 


37 


96 


;o 


15) 


E 


1992br 


4.418 


0.40 


0.00 


38 


21 


;o 


19) 


1.77 


0.02 


38 


09 


;o 


36) 


E 


1992bs 


4.283 


0.00 


0.00 


37 


61 


;o 


14) 


1.10 


0.24 


37 


63 


;o 


18) 


L 


1993H 


3.871 


0.16 


0.67 


35 


20 


;o 


26) 


1.59 


0.72 


35 


23 


;o 


25) 


L 


19930 


4.189 


0.03 


0.00 


37 


03 


;o 


12) 


1.18 


0.00 


37 


31 


;o 


14) 


E 


1993ag 


4.177 


-0.19 


0.64 


36 


80 


;o 


17) 


1.29 


0.55 


37 


11 


;o 


19) 


E 



* Light curves resticted to B and V data within 40 days of B maximum 
*E=early-type host, L=late-type host 



